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I.  INTRODUCTION  AND  SUMMARY 


1. 1 Introduction 

Commonly  used  methods  for  estimating  the  yield  of 
foreign  underground  nuclear  explosions  rely  on  the  use  of 
empirical  relations  between  body  (m^)  and  surface  (M  ) wave 
magnitudes  and  yield.  The  difficulty  with  empirical  relation- 
ships is  that  they  can  be  systematically  in  error  when  applied 
to  events  outside  the  empirical  data  base.  This  problem  is 
widely  recognized  and,  as  a result,  there  is  considerable 
interest  in  developing  a clear  understanding  of  the  important 
parameters  that  control  the  seismic  signatures  of  underground 
explosions . 

In  this  report  we  describe  a study  of  NTS  explosions 
in  which  we  have  attempted  to  model  in  great  detail  the 
recorded  Rayleigh  waves  at  two  WWSSN  stations:  Albuquerque, 
New  Mexico  (ALQ)  and  Tucson,  Arizona  (TUC) . We  find  that 
our  synthetic  waveforms  are  in  excellent  agreement  with  the 
observations.  This  suggests  that  most  of  the  important 
contributing  factors  are  properly  represented  by  our  models. 

The  main  objectives  of  the  work  described  in  this 
report  are  to  use  the  comparison  of  synthetic  and  observed 
Rayleigh  waves  to  infer  the  long  period  amplitude  of  the  NTS 
explosion  source  and  to  explore  the  potential  influence  of 
supposed  source  complexities,  particularly  surface  spallation 
and  related  phenomena.  In  pursuing  these  objectives  we 
address  the  following  questions: 

•'How  accurately  can  the  Rayleigh  wave  signature  of 
NTS  explosions  be  modeled  with  plane-layered  elastic 
models  for  the  travel  path? 

• Is  a spherically  symmetric  point  source  adequate 
for  modeling  the  source? 


What  is  the  effect  of  surface  spallation  on  the 
Rayleigh  wave? 

• What  is  the  effect  of  tectonic  strain  release  on 
the  Rayleigh  wave? 

• What  is  the  source  amplitude  and  how  does  it  depend 
on  source  material,  source  depth  and  yield? 

• How  much  scatter  is  to  be  expected  for  events  in 
relatively  homogeneous  regions  (e.g.,  below  the 
water  table  at  Yucca  Flat)  and  what  are  the  main 
causes  of  this  scatter? 

We  have  chosen  to  use  data  from  the  stations  ALQ  and 
TUC  because  the  data  was  conveniently  available  and  because 
these  stations  are  at  an  excellent  range  (700  to  950  km)  for 
studying  NTS  explosions.  They  are  close  enough  to  record 
many  of  the  smaller  yield  events,  while  being  at  a range 
where  the  trace  is  dominated  by  the  fundamental  mode  Rayleigh 
wave.  At  their  normal  gain  setting  the  Rayleigh  wave  trace 
can  be  measured  for  yields  as  large  as  200  to  300  KT. 

The  first  step  in  our  analysis  of  the  ALQ  and  TUC 
data  was  to  determine  the  crustal  structure  along  the  two 
paths  and  this  work  was  reported  by  Bache,  Rodi  and  Harkrider 
(1978).  The  crustal  models  were  determined  by  inverting  the 
phase  and  group  velocities  inferred  directly  from  the  ALQ  and 
TUC  recordings  of  NTS  explosions.  The  inverted  crustal  struc- 
tures were  satisfying  in  that  they  were  simple  (few  major 
discontinuities  were  present) , agreed  with  other  available 
data  (e.g.,  from  refraction  studies)  for  these  paths  and 
gave  an  excellent  fit  to  the  dispersion  data.  Using  attenuation 
models  based  on  data  given  by  Mitchell  (1975),  Bache,  Rodi 
and  Harkrider  (1978)  computed  synthetic  seismograms  for  repre- 
sentative events  in  several  test  areas  at  NTS  and  found  that 
the  agreement  of  theoretical  and  observed  waveforms  was 
very  good. 
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1.2 


Outline  of 


the  Analysis 


The  results  of  the  earlier  work  gave  reason  to  believe 
that  we  could  quite  accurately  account  for  the  path  effects 
for  NTS  explosions  observed  at  ALQ  and  TUC  and  thus  determine 
the  characteristics  of  the  source.  This  is  the  primary 
objective  of  the  work  described  in  this  report.  In  outline 
form  our  report  includes  the  following: 


A.  The  pertinent  results  of  the  Bache,  Rodi  and 
Harkrider  study  are  summarized  (Section  II) . 

B.  The  theoretical  formulation  used  to  compute 
synthetic  Rayleigh  wave  seismograms  is  summarized 
with  particular  attention  to  the  parameters  that 
control  the  amplitude  (Section  III)  . 

C.  A spherically  symmetric  source  represented  by  a 
reduced  displacement  potential  is  assumed.  Then 
for  each  of  the  test  areas  studied  (Yucca  Flat, 
Pahute  Mesa  and  PILEDRIVER) , the  Rayleigh  wave 
amplitude  is  directly  proportional  to  the  static 
value  of  the  reduced  displacement  potential  (¥oo)  . 
Comparing  theoretical  and  observed  seismograms , 

a is  determined  from  each  station  for  each 
event  and  the  values  are  scaled  to  a common  yield. 
The  important  results  are  as  follows  (Section  IV) : 

• At  ALQ  the  mean  values  of  the  ^ scaled  to 
0.02  KT  are  9 . 1 m^  for  the  17  Yucca  Flat  events 
6.6  ra-*  for  the  7 Pahute  Mesa  events  and  3.9  for 
PILEDRIVER.*  Only  events  below  the  water  table 
with  yields  from  40  to  200  KT  were  studied. 

The  standard  deviation  for  each  population  is 
about  40  percent  of  the  mean. 

• The  4'0O  inferred  from  the  TUC  records  are  con- 
sistently 1.5  times  larger  than  from  ALQ. 


# 

Throughout  the  report  we  use  the  notation  ¥ for  both  the 
absolute  value  for  a given  event  and  for  values  scaled  by 
the  cube-root  of  the  yield  to  0.02  KT.  When  there  is  any 
possibility  of  confusion  about  which  kind  of  value  is  meant, 
we  use  the  notation  ¥-02  to  indicate  the  common  yield  value. 


f 


• For  similar  materials  the  Rayleigh  waves  are 
proportional  to  us  where  us  is  the  average 
shear  modulus  in  the  source  region.  Variations 
in  us  and  errors  in  the  official  yield  are  two 
sources  of  random  scatter  in  the  ¥ in  a 
particular  area. 

• Random  effects  can  plausibly  account  for  at 
most  half  the  scatter  in  the  inferred  source 
levels  for  each  population.  The  remainder  of 
the  scatter  appears  to  be  due  to  real  differences 
among  the  source  levels  of  the  different  events 
in  each  population. 

• The  slope  of  the  log  y^-log  yield  curve  (com- 
parable to  Ms-log  yield)  cannot  be  determined 
with  confidence  for  the  rather  narrow  yield 
range  represented  in  our  data. 

The  4^  values  from  this  analysis  of  Rayleigh  waves 
are  compared  to  estimates  of  4'00  made  by  other 
methods;  for  example,  from  close-in  methods,  from 
far-field  body  and  surface  waves  or  from  theoretical 
considerations.  Our  values  are  within  the  range 
expected  from  this  other  work  except  for  PILEDRIVER 
where  our  value  is  low  (Section  IV) . 

The  assumption  that  the  source  is  spherically 
symmetric  is  clearly  an  oversimplification.  Using 
results  obtained  by  Toksoz  and  Kehrer  (1972),  we 
correct  our  solution  for  the  effect  of  a double- 
couple component  in  the  source.  We  find  that  the 
double-couple  has  virtually  no  effect  on  the  wave- 
form, but  simply  scales  the  amplitude.  If  the 
values  given  by  Toksoz  and  Kehrer  are  typical, 
some  15  to  20  percent  of  the  discrepancy  between 
the  ALQ  and  TUC  solutions  for  the  Yucca  Flat  and 
Pahute  Mesa  4^  is  due  to  the  double-couple.  The 
mean  4^2  values  at  ALQ  are  reduced  from  9.1  and 
6.6  to  8.1  and  4.3  for  these  two  areas.  For 
PILEDRIVER  Toksoz  and  Kehrer  predict  a much  larger 
effect  on  the  Rayleigh  waves  at  these  two  stations. 
Their  solution  actually  increases  the  discrepancy 
between  the  values  obtained  at  the  two  stations, 
though  it  is  quite  sensitive  to  small  errors  in 
the  double-couple  orientation.  Further,  the 
double-couple  contribution  so  dominates  the  solu- 
tion that  the  required  to  match  the  data  at 

ALQ  is  reduced  to  about  1.0  (Section  V). 
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F.  The  impulse  delivered  by  the  impact  of  the  large 
mass  of  spalled  material  returning  to  the  free 
surface  has  been  supposed  to  have  some  significant 
effect  on  the  Rayleigh  wave  signature  of  explosions 
(e.g„,  Viecelli,  1973).  Using  estimates  for  the 
spall  impulse  given  by  Viecelli  (1973)  and  somewhat 
larger  estimates  given  by  Sobel  (1978),  we  compute 
expected  Rayleigh  waves  at  ALQ  and  TUC  to  determine 
its  effect.  We  find  it  to  be  quite  small.  It 
seems  implausible  to  suppose  that  the  spall  impulse 
could  change  the  Rayleigh  wave  amplitude  by  more 
than  5 to  10  percent  (Section  VI) . 

G.  A spall-related  phenomenon  that  may  be  much  more 
significant  is  the  associated  loss  of  energy 
from  the  waves  traveling  upward  from  the  source. 

We  explore  the  potential  effect  by  simply  deleting 
a portion  of  the  upgoing  waves  from  the  solution. 

A more  appropriate  filter  may  be  frequency-dependent, 
having  its  primary  effect  on  the  short  periods, 
but  at  present  we  have  no  results  to  support  this 
conjecture.  The  events  at  Yucca  Flat  are  quite 
insensitive  to  the  suppression  of  the  upgoing 
waves.  The  effect  at  Pahute  Mesa  is  much  larger 
with  suppression  of  half  the  upgoing  waves  leading 
to  a decrease  of  about  25  percent  in  the  amplitude. 
PILEDRIVER  is  extremely  sensitive  to  the  amount 
of  upgoing  waves  suppressed  and  this  is  a very 
important  parameter  for  events  in  granite.  If  we 
suppose  that  50  percent  of  the  upgoing  wave  is 
lost  to  spallation  or  scattering,  the  Rayleigh 
wave  amplitude  is  reduced  by  about  40  percent. 
Suppression  of  25  percent  of  the  upgoing  waves 
caused  a decrease  of  about  20  percent  in  the 
amplitude  (Section  VI) . 

H.  We  carefully  examine  the  factors  that  control 
the  amplitude  of  the  synthetic  seismograms  to 
determine  why  the  source  level  inferred  from  the 
TUC  observations  is  about  1.5  times  larger  than 
that  inferred  from  the  ALQ  observations.  A partial 
explanation  is  provided  by  the  presence  of  asymmetries 
at  the  source,  particularly  the  double-couple  (E.), 
but  a substantial  amount  remains.  There  are  four 
main  features  of  the  theoretical  calculation  that 

are  examined  closely?  the  dispersion,  the  model 
for  the  source  region,  the  transition  between  the 
local  source  structure  and  average  path  structure 
and  the  inelastic  attenuation  (Q)  for  the  two 
paths.  Rather  than  from  any  one  factor,  the 
differences  in  the  source  estimates  seem  to  be 
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due  to  contributions  from  each  of  these.  Errors 
in  the  dispersion  and  attenuation  are  small,  but 
tend  in  the  riqht  direction.  The  main  source  of 
error  is  probably  associated  with  the  failure  of 
plane-layered  laterally  homogeneous  models  to  pre- 
cisely represent  the  complex  real  earth.  We  recom- 
mend that  the  answers  from  the  two  stations  be 
averaged,  but  give  more  weight  to  those  from  ALQ. 

A better  definition  of  the  true  source  amplitude 
would  probably  be  achieved  by  carrying  out  the  same 
analysis  for  more  stations  and  averaging  the  re- 
sults (Section  VII). 


v 


1. 3 Summary  of  Results 

In  our  analysis  we  divided  the  explosions  into  three 
groups  (Yucca  Fiat,  Pahute  Mesa  and  PILEDRIVER)  separated 
geographically  and  by  the  average  material  properties  at  the 
source.  Using  plane-layered  earth  models,  synthetic  seismograms 
were  computed  that  give  excellent  agreement  with  the  observed 
waveforms.  For  the  source  we  used  four  different  models  that 
were  combined  in  various  ways.  These  are: 

• A spherically  symmetric  point  source  given  by  a 
reduced  displacement  potential. 

• A reduced  displacement  potential  source  with  a 
portion  of  the  upgoing  waves  suppressed. 

• A downward  impulse  representing  the  impact  of  the 
spalled  material. 

• A double-couple. 

The  source  quantity  most  directly  related  to  explosion 
yield  is  the  static  level  of  the  reduced  displacement  potential 
0?  ) .*  We  first  estimate  this  quantity  assuming  the  source  is 
spherically  symmetric,  then  correct  these  estimates  for  other 
effects.  Our  best  estimates  for  the  mean  scaled  to  0.02 
KT  in  each  area  are  as  follows: 

♦ The  relationship  between  'fa  and  yield  is  dependent  on  the 
local  material  properties  and  is  a subject  for  separate  dis- 
cussion (e.g.,  Cherry,  et  al. , 1975) 
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Yucca  Flat  Pahute  Mesa  PILEDRIVER 


1.  Assume  spherical 
symmetry  (Table  4) 

2.  Correct  for 
observed  double- 
couple (Table  8) 

3.  Correct  for 
spallation* 
(Figures  21-23) 

4.  Simultaneously 
correct  for 
spallation  and 
double-couple 
(Section  6.5) 


10.5 


8.3 


11.4 


9.1 


8.1 


4.8 


10.0 


5.8 


4.7 

1.4 

5.9 


1.8 

2.  7 


* The  spall  impact  makes  little  contribution.  The  main  effect 
is  the  suppression  of  upgoing  waves  from  the  source.  This 
estimate  is  based  on  suppression  of  75  percent  of  these  waves. 

t 25  percent  of  the  upgoing  waves  from  both  the  explosion  and 
double-couple  are  suppressed.  For  Yucca  Flat  and  Pahute 
Mesa,  the  suppression  is  only  applied  to  the  explosion  waves. 


Throughout  the  report  values  are  determined  separately  for  the 
two  stations  ALQ  and  TUC.  The  "best"  estimates  given  above 
were  determined  by  averaging  the  two  with  the  ALQ  values 
given  double  weight. 

The  values  given  represent  mean  values  and  individual 
events  can  deviate  substantially  from  these.  The  corrections 
for  the  double-couple  and  spallation  contributions  are  again 
based  on  mean  values,  but  these  effects  must  vary  widely  from 
event  to  event.  In  fact,  variation  in  these  effects  is  likely 
to  be  responsible  for  a large  portion  of  what  we  have  called 
"real"  source  level  variations;  that  is,  variations  that  cannot 
be  attributed  to  random  errors  in  our  modeling  procedure. 

The  effect  of  the  double-couple  can  be  reduced  by  averaging 
values  from  many  azimuths,  while  the  spallation  effects  are  not 
likely  to  depend  on  azimuth. 
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This  report  is  concerned  with  mean  source  amplitudes 
for  selected  populations  of  events.  The  values  for  individual 
events  are  available  for  those  with  the  proper  security 
clearance  (Bache,  1978) . A substantial  amount  of  the  deviation 
of  individual  events  from  the  mean  may  be  explicable  from 
known  peculiarities  of  these  events.  We  have  not  attempted 
to  correlate  our  values  with  candidate  parameters  like  the 
velocities  or  other  material  properties  near  the  event.  This 
may  be  a useful  exercise  for  the  future. 


1. 4 Conclusions 

Our  major  conclusions  are  summarized  in  Section  VIII. 
Most  important  is  our  demonstration  that  available  techniques 
are  capable  of  modeling  surface  wave  signatures  of  underground 
explosions  in  considerable  detail.  The  amplitude  of  the  source 
can  thus  be  inferred  from  a comparison  of  synthetic  and  ob- 
served seismograms.  Applied  to  foreign  explosions,  this  would 
lead  to  yield  estimates  that  are  independent  of  the  estimates 
from  empirical  Ms~yield  curves. 
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II.  COMPARISON  OF  THEORETICAL  AND  OBSERVED  RAYLEIGH 

WAVEFORMS 

Typical  ALQ  and  TUC  recordings  of  NTS  explosions  are 
shown  in  Figure  1,  and  the  pertinent  data  for  these  events 
are  summarized  in  Table  1 (Springer  and  Kinnaman,  1971) . 

There  is  a clearly  defined  maximum  phase  on  each  record. 

Within  the  test  area  groupings  (PILEDRIVER,  Pahute  Mesa, 

Yucca  Flat)  the  waveforms  are  remarkably  consistent  for  the 
more  than  50  events  examined  at  each  station.  This  makes  it 
possible  to  select  representative  events  from  each  group  for  j 

closer  analysis.  f 

i 

Bache,  Rodi  and  Harkrider  (1978)  inferred  crustal 
models  for  the  NTS-TUC  and  NTS -ALQ  paths.  This  was  done  in 
the  following  way: 

1.  Seismograms  for  PILEDRIVER,  DURYEA  and  TAN  at 
each  station  (Figure  1)  were  hand-digitized. 

2.  The  phase  and  group  velocities  were  derived 
from  the  seismograms.  The  results  are  shown 
in  Figure  2. 

3.  Using  linear  inversion,  crustal  structure  were  i 

inferred  from  the  phase  and  group  velocity  data.  j 

The  crustal  structures  were  also  constrained  to 

be  consistent  with  other  information  about 
the  paths.  The  structures  and  the  agreement 
observed  and  theoretical  dispersion  are  shown 
in  Figures  3 and  4. 

4.  Using  methods  outlined  in  the  next  section, 
theoretical  seismograms  were  computed  for  each 
station-event  pair.  The  comparison  to  the  ob- 
servations is  shown  in  Figure  5. 
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Density  p (gm/cm3) 


Trt  s249  s 


ViTl 


PILEDRIVER  (Climax  Slock) 


T0  = 203  s 


DURYEA  (Pahute  Mesa) 


T0=234  s 


HI 


TAN  (Yucca  Flat) 


Mi  nil 

Hill'll 


one  minute 


Figure  5 . 


Theoretical  and  observed  seismograms  are  compared 
at  ALQ  (left)  and  TUC  for  events  in  three  test 
areas  at  NTS.  A bar  indicating  one  minute  is 
shown.  In  each  pair  the  observed  (top)  and  theo- 
retical records  start  at  the  same  time  with  respect 
to  the  explosion  detonation  and  this  time  is  indi- 
cated as  T . 
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The  theoretical  and  observed  seismograms  agree  very 
well  at  the  two  stations,  even  in  some  rather  subtle  details. 
Agreement  of  theoretical  and  observed  seismograms  requires 
not  only  that  the  dispersion  be  matched,  but  that  the  amplifi- 
cation due  to  the  source  and  travel  path  be  properly  repro- 
duced. 


The  Airy  phase 

amplitudes 

and  periods 

can  be  measured 

very  accurately  on  the 

synthetic 

seismograms . 

The  periods 

are  as  follows: 

ALQ 

TUC 

PILEDRIVER 

10.4 

8.4 

DURYEA 

10.7 

7.0 

TAN 

9.9 

7.6 

These  are  close  to  the  periods  measured  (much  less  accurately) 
on  the  observed  records  which  are  11  ±0.5  seconds  for  all 
events  at  ALQ  and  8 ±0.5  seconds  at  TUC.  The  instrument 
amplification  is  0.78  at  7.0  seconds,  0.86  at  8.4  seconds 
and  0.95  at  10.7  seconds.  Therefore,  errors  in  the  period 
measurement  do  not  have  much  effect  on  the  amplitude  deter- 
mination. In  later  sections  the  amplitude  comparison  will 
be  used  to  infer  the  source  level. 
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III.  PARAMETRIC  DEPENDENCE  OF  RAYLEIGH  WAVE  AMPLITUDES 


The  theoretical  relationships  used  in  computing  the 
synthetic  seismograms  of  Figure  5 were  given  by  Bache,  Rodi 
and  Harkrider  (1978) . The  basic  formula  for  the  vertical 
component  of  the  fundamental  mode  Rayleigh  wave  at  range  r from 
a spherically  symmetric  explosion  is 


K 


w(r,w)  = - 4it 


ligHw) 


S (h)^R 
S1  R1 


TU)  H 


(2) 


(if) 


,-Yr 


1/2 


;inA  j 


(1) 


where  all  quantities  are  Fourier  transformed  with  respect  to 

time.  The  u is  the  shear  modulus  in  the  source  region,  c is 
s ( 2 ) 

Rayleigh  phase  velocity  and  is  a Hankel  function.  The 

Ar  is  the  source  depth-independent  excitation  factor  and 

K (h)  is  the  depth-dependent  excitation.  The  quantity 
S1 


>*Yr 


a sinA 
e 


1/2 


accounts  for  anelastic  attenuation  and  the  sphericity  of  the 
earth . 

A two  structure  model  is  used  for  the  path  with  sub- 
script 1 in  (1)  denoting  the  source  structure  and  2 denoting 
what  we  call  the  "path"  structure.  Note  that  the  dispersion 
is  entirely  controlled  by  the  path  structure  and  the  ampli- 
tude excitation  is  controlled  by  the  source  structure. 
Passage  of  Rayleigh  waves  between  the  two  structures  is  ac- 
counted for  by  the  transmission  coefficient  T(oo). 

Finally,  we  have  the  source  which  is  represented  by 

A 

'F  ( to ) in  (1).  For  spherically  symmetric  explosions  this  is 
related  to  the  displacement  by 
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u (R,oj) 


l (hi)  . $ (u) 
R2  Ra 


e~iajR/a 


(2) 


where  R is  the  range  and  a is  P wave  velocity. 

We  now  examine  the  various  factors  that  control  the 
amplitude  of  the  theoretical  seismograms  in  Figure  5.  The 
source  structure  was  constructed  for  each  of  the  three 
examples  by  altering  the  top  two  kilometers  of  the  NTS-TUC 
path  model  to  represent  Pahute  Mesa,  Yucca  Flat  or  the 
PILEDRIVER  site.  The  velocity-density  profiles  for  this 
portion  of  the  models  are  shown  in  Figure  6. 

From  (1),  the  Rayleigh  waves  from  explosions  are 
proportional  to 


I ! 


u 


s 


(3) 


For  the  periods  (>  4 seconds)  and  yields  (50-500  KT)  of 

A 

interest,  'M x)  should  be  nearly  equal  to  a constant  which 
we  call  Our  primary  objective  here  is  to  deduce  the 

value  of  this  V . 

00 

The  source  excitation,  u Ks  ar ,/c, , and  the  trans- 
mission  coefficient,  TU),  are  plotted  in  Figure  7. 

We  note  that  T(ui)  is  near  unity  for  the  periods  (z  8 seconds 
for  TUC  and  : 11  seconds  for  ALQ,  see  Section  II)  of  primary 
interest. 

The  relative  source  excitation  terms  for  Pahute  Mesa 

and  Yucca  Flat  differ  by  a nearly  constant  factor  over  most 

of  the  frequency  range.  Some  of  these  values  and  their 

ratio  are  given  in  Table  2 . The  ratio  is  little  different 

from  the  ratio  of  the  shear  moduli  (u  ) which  is  3.07.  This 

s 

is  important  because  it  says  that  the  Rayleigh  waves  from 
explosions  in  the  tuffs  and  rhyolites  at  Yucca  Flat  and 

A 

Pahute  Mesa  are  very  nearly  proportional  to  u ^(w). 

s 
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Figure  7.  The  source  excitation  factor  is  shown  at  left  for  the  three  s 
areas  studied.  At  right  is  the  transmission  coefficient  T(oo) 
the  six  source-path  combinations. 


TABLE  2 


VALUES  OF  THE  EXCITATION  FACTOR  u K A„  /c 


Yucca  Flat 

Period  (s ec/cm?) 


16 

8.11 

X 

10-14 

14 

1.01 

X 

IQ'13 

12 

1.27 

X 

io-13 

10 

1.68 

X 

10-13 

8 

2.45 

X 

io-13 

6 

4.30 

X 

io"13 

1 1 

Pahute  Mesa 


(sec/cm2) 

Ratio 

-13 

2.53  x 10  J 

3.13 

3.12  x 10*13 

3.10 

3.89  x 10*13 

3.06 

5.05  x 10~13 

3.00 

7.14  x IO"13 

2.91 

12.2  x 10~L3 

2.83 

21 
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The  PILEDRIVER  site,  Climax  Stock  granite,  is  suf- 
ficiently different  from  the  other  two  that  the  proportionality 
to  ug  is  not  applicable.  In  fact,  the  amplitude  factor  is 
much  smaller  than  expected  from  the  large  ug  at  this  site. 

The  remaining  factors  in  (1)  are  the  Hankel  function, 
which  accounts  for  the  dispersion  and  elastic  attenuation, 
the  anelastic  attenuation  factor  e Yr  and  the  sphericity 
correction.  For  the  ranges  of  interest,  the  latter  is  near 
unity.  The  dispersion  has  a strong  effect  on  the  maximum 
amplitude  because  it  controls  the  interference  between  the 
frequency  components. 

Our  attenuation  model  for  the  calculations  of  Figure 
5 is  based  on  Western  U.  S.  data  compiled  by  Mitchell  (1975) . 
Since  these  data  were  derived  from  observations  of  events  on 
the  Colorado  Plateau,  they  were  assumed  to  be  appropriate  for 
the  NTS-ALQ  path.  Using  Mitchell's  y values,  we  constructed 
Q models  for  the  two  paths,  assuming  the  Q-$  relationship 
was  the  same  for  both.  These  Q models  are  shown  in  Figure  8. 
The  Y(i>)  which  result  from  these  models  are  tabulated  in 
Table  3.  Also  shown  in  Figure  8 are  the  e_Yr  for  the  two 
paths,  using  representative  distances. 

In  summary,  the  factors  controlling  the  maximum  peax- 
to-peak  amplitude  of  the  theoretical  Rayleigh  waves  at  ALQ 
and  TUC  are: 

1.  The  dispersion  characteristics  of  the  path 

model  — from  Figure  3 we  see  that  11  seconds, 
the  period  of  the  maximum  phase  on  the  ALQ  records, 
occurs  where  the  group  velocity  is  slowly  de- 
creasing. At  8 seconds  on  the  TUC  dispersion, 
Figure  4,  the  group  velocity  is  near  a local 
maximum.  The  exact  shape  of  the  dispersion 
curves  near  the  dominant  period  will  influence 
the  amplitude. 
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TABLE  3 


1 


ATTENUATION  PARAMETERS,  y (m)  , 
SEISMOGRAMS  OF 

FOR  THE  THEORETICAL 
FIGURE  5 

Period 

Y,  NTS-ALQ 

Y,  NTS-TUC 

25 

1.58  (x  10"4  km"1) 

1.75  (x  10-4  km' 

22.5 

1.74 

1.89 

20 

1.96 

2.05 

18 

2.24 

2.22 

16 

2.67 

2.47 

14 

3.33 

2.89 

12 

4.34 

3.59 

10 

5.94 

4.80 

9 

7.12 

5.71 

8 

8.70 

6.96 

7 

11.0 

8.74 

6 

14.5 

11.4 

5 

21.4 

15.9 

4 

43.5 

24.2 

3 

183.0 

43.1 

2 

401.0 

114.0 

I 


I 
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2.  The  structure  in  the  source  region  — this  controls 
the  source  excitation  factors  (h)  , Ar^  and  c^. 
The  shear  velocity  at  the  source  is  also  a key 
parameter . 


3.  The  transmission  coefficient,  T(w)  — this  quantity 
accounts  for  the  passage  of  the  Rayleigh  wave  from 
the  source  structure  into  the  path  structure. 


4.  The  parameter  y(u)  that  accounts  for  anelastic 
attenuation  and  elastic  scattering. 


5.  The  source  function  --  Equation  (1)  assumes  the 

source  is  a spherically  symmetric  explosion.  How- 
ever, the  source  may  include  other  components. 

A strong  contribution  to  the  Rayleigh  waves  may 
be  present  due  to  some  form  of  tectonic  release 
(Toksoz  and  Kehrer,  1972;  Lambert,  et:  al . , 1972) 
or  spall  slapdown  (Viecelli,  1973). 


Our  purpose  is  to  discover  the  amplitude  and  nature  of 
the  source  function.  In  later  sections  we  will  attempt  to 
bound  the  contribution  from  the  other  four  factors  listed 
above . 
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IV.  SOURCE  ESTIMATES  ASSUMING  SPHERICAL  SYMMETRY 

Assuming  the  source  is  a spherically  symmetric  explosion 
and  the  path  is  modeled  by  the  plane-layered  earth  models  de- 
scribed in  the  previous  section,  we  computed  the  theoretical 
seismograms  shown  in  Figure  5.  Since  the  other  factors  are 
fixed,  the  amplitude  of  the  theoretical  seismograms  can  be  j 

scaled  by  the  source  function,  f,  in  Equation  (1).  The 
theoretical  seismograms  shown  in  Figure  5 were  computed  at  a 
depth  of  500  meters  and  a yield  of  100  KT  using  the  source 
function  shown  in  Figure  9.  This  is  source  133  from  Bache, 
et  a_l.  (1975)  with  the  frequency  axis  scaled  to  100  KT.  This 
source  function  was  computed  for  saturated  Yucca  Flat  tuff, 
but  its  shape  is  representative  of  that  deduced  in  numerous 
theoretical  and  empirical  studies  of  NTS  explosions. 

The  amplitude  axis  in  Figure  9 is  scaled  to  0.02  KT. 

For  any  other  yield  (W)  the  amplitude  is  obtained  by  multiplying 
by  W/0.02.  Note  that  the  source  function  is  nearly  constant 
over  the  frequency  band  of  interest.  The  source  level  for  5 
second  waves  is  about  15  percent  larger  than  for  very  long 
period  waves.  In  view  of  this  small  deviation  from  a constant 
level,  it  is  convenient  to  characterize  the  source  by  f , the 

T 00 

zero  frequency  limit,  which  is  10.3  m for  0.02  KT. 

In  Figures  10  and  11,  time-domain  Airy-phase  amplitudes 
observed  at  ALQ  and  TUC  are  plotted  versus  yield  for  repre- 
sentative events  with  yields  near  100  KT.  On  each  plot, 

• 02 
theoretical  unit-slope  lines  indicating  constant  Hr00  are 

also  shown.  The  values  of  4^^  labeling  these  lines  were 

determined  from  the  synthetic  seismograms  in  Figure  5.  That 

.02 

is,  if  the  for  the  synthetic  seismogram  calculations  had 
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Figure  9.  The  amplitude  of  the  source  function  for  the  theoretical 
seismograms  is  plotted  with  the  frequency  axis  scaled 
to  100  KT  and  the  amplitude  axis  to  0.02  KT.  The 
amplitude  values  at  several  periods  are  indicated  on 
the  plot. 
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these  values,  then  their  Airy-phase  amplitudes  as  a function 
of  yield  would  be  described  by  the  lines  shown.  Neglecting 
the  small  correction  for  the  yield  scaling  of  the  frequency 

/■v 

dependence  of  f(oj),  these  lines  have  unit  slope. 

The  lines  in  Figures  10  and  11  can  be  used  to  deduce 
the  (0.02  KT)  source  level  required  to  bring  theoretical  and 
observed  Airy  phase  amplitudes  into  coincidence.  Note  that, 
on  the  average,  a larger  source  function  seems  to  be  required 
to  fit  the  TUC  data  than  to  fit  the  ALQ  data.  We  will  later 
discuss  this  apparent  asymmetry  at  some  length.  However,  we 
first  develop  a more  quantitative  description  of  the  inferred 
source  amplitudes. 

Within  localized  source  regions  like  Pahute  Mesa  and 
Yucca  Flat,  there  are  several  effects  that  might  cause  the 
ratio  of  Airy  phase  amplitude  to  explosion  yield  to  vary  from 
event  to  event  and  thus  account  for  the  scatter  observed  in 
Figures  10  and  11.  These  include: 

1.  Range  variations, 

2.  Depth  of  burial  variations, 

3.  Material  rigidity,  us,  variations 

4.  Variations  in  the  material  properties  leading  to 

variations  in  ¥ /KT. 

00 

5.  Variations  in  the  frequency  dependence  of  the 
source  function,  primarily  the  ratio  of  the 
source  level  at  the  Airy  phase  period,  to  that 
at  long  period. 

These  items  are  based  on  the  assumption  that  the  source  is 
spherically  symmetric.  Asymmetric  effects  may  also  cause 
variations,  but  we  will  discuss  these  later. 
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We  now  examine  the  factors  listed  above  to  see  what 
' their  effects  can  be.  First,  for  Yucca  Flat  events,  the 

distance  to  TUC  varies  from  708  to  724  km,  while  for  Pahute 
Mesa  events,  the  TUC  range  variation  is  747  to  764  km.  To 
;j  quantify  this  effect,  synthetic  seismograms  were  computed 

| for  the  Yucca  Flat  event  TAN  at  ranges  of  708  and  724  km. 

The  amplitude  variation  was  about  8 percent.  Similar  cal- 
ij  culations  were  done  for  the  Pahute  Mesa  event  DURYEA  at  74  7 

and  763.5  km.  In  this  case  the  Airy  phase  amplitude  variation 

was  less  than  four  percent.  In  Figure  12  the  seismograms  at  the  j 

J new  distances  are  compared  to  the  TAN  and  DURYEA  synthetics  ■ 

? 

[;  from  Figure  5.  We  see  that  the  waveforms  vary  little  with 

. small  changes  in  range. 

i The  amplitude  variation  with  depth  was  determined  by 

! comparing  theoretical  seismograms  computed  with  only  the  depth 

j changing.  For  the  depths  of  interest  we  find  the  following: 

| Yucca  Flat:  log  A **  0.16  log  H,  ^ 

| Pahute  Mesa:  log  A =&  0.05  log  H, 

where  H is  the  depth  of  burial.  If  H * w1//3,  the  yield  scaling 
becomes  log  A = 1.0  5 log  W for  Yucca  Flat  and  log  A **  1.0  2 
log  W for  Pahute  Mesa. 

The  material  rigidity  (y  ) variation  can  be  quite  im- 
portant. As  was  pointed  out  in  Section  III,  the  Rayleigh  wave 
amplitude  is  directly  proportional  to  y . Of  course,  shear 
modulus  variations  are  likely  to  be  associated  with  changes 
in  the  coupling  characteristics  of  the  material;  that  is,  in 

'i 

the  ¥(0}).  Such  material  properties  as  the  water  content,  air- 
filled  porosity  and  material  strength  influence  the  seismic 
coupling  (Cherry,  Rimer  and  Wray,  1975) . 

We  should  also  point  out  that  the  absolute  level  of 
the  ^ we  infer  is  especially  sensitive  to  the  yg  used  in  our 
synthetic  seismogram  calculations  (see  Figure  7).  Without 
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Range  (km)  h Range  (km) 


independent  knowledge  of  the  shear  modulus,  we  can  only 
determine  us  'Ha>.  In  subsequent  discussion  we  will  continually 
refer  to  the  inferred  ¥ values,  but  the  trade-off  with  our 
assumed  u should  be  kept  in  mind. 

The  frequency  dependence  of  the  source  function,  item 
5 in  the  above  list,  is  difficult  to  quantify.  Generally, 
we  would  expect  the  source  function  to  be  nearly  flat  for  low 
yield  explosions.  However,  for  very  large  yields  (hundreds 
of  kilotons)  the  source  level  at  8 to  10  seconds  could  be 
greater  than  at  larger  periods. 

Estimates  of  ¥ were  made  for  each  event  in  the  fol- 
00 

lowing  way.  The  amplitudes  of  the  synthetic  seismograms  of 
Figure  5 were  corrected  for  source  depth  and  yield  by  multi- 
plying by 


where  the  n is  0.16  or  0.05  from  (4).  The  small  correction 
for  range  variations  was  not  made.  The  event  values, 
scaled  to  0.02  KT,  are  then  the  from  Figure  9 (10.3)  times 
the  ratio  of  the  observed  amplitude  to  the  corrected  synthetic 
amplitude.  The  theoretical  source  level  is  a bit  higher  at 
the  8 (TUC)  and  11  second  (ALQ)  periods  at  which  the  measurements 
are  made,  but  as  long  as  we  realize  that  the  numbers  are 
based  on  the  particular  source  function  shown  in  Figure  9, 
this  should  cause  no  problem. 

The  were  computed  from  the  observations  at  ALQ 
and/or  TUC  for  the  following  events: 

Yucca  Flat:  PIRANHA,  BRONZE,  GRAPE  A,  SHAPER,  TAN 

MINATA,  STARWORT,  CUP,  CALABASH,  OSCURO, 

GRAPE  B,  LAMPHER,  NOGGIN,  DUMONT,  BUFF, 

MIERA,  AGILE,  CORDUROY,  TIJERAS,  ZAZA 
and  KNOX; 
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Pahute  Mesa:  DURYEA,  KNICKERBOCKER,  CHATEAUGAY, 

SCOTCH,  STINGER,  PURSE,  SLED; 

PILEDRIVER. 

The  inferred  values  are  plotted  versus  explosion 

yield  in  log-log  form  in  Figure  13.  For  each  population  the 
logarithmic  mean  values  and  their  95  percent  confidence  limits 
were  computed  using  Student's  t distribution.  These  values 
are  shown  on  the  plot  and  are  summarized  in  Table  4.  If  the 
source  were  spherically,  or  at  least  axially,  symmetric  and 
our  theory  properly  accounted  for  the  path,  the  estimates 
would  be  the  same  at  ALQ  and  TUC  for  each  event.  However,  the 
source  level  required  to  match  the  TUC  amplitudes  appears  to  be 
significantly  larger  than  that  required  to  match  the  ALQ  amplitudes 

Since  our  data  span  only  a narrow  yield  range,  we  are 
unable  to  determine  the  dependence  of  the  ^ values  on  yield 
(or  depth  of  burial,  since  the  two  are  nearly  proportional) . 

This  is  demonstrated  in  Table  5 by  the  results  of  a linear 
regression  fit  to  the  data  plotted  in  Figure  13.  If  cube-root 
scaling  applies,  the  true  is  proportional  to  W and  the  slope 
of  02  versus  W would  be  zero.  For  all  four  sets  the  slope 
is  negative  suggesting  that  W*  with  n < 1.  However,  the 

confidence  limits  are  quite  wide  and  we  can  have  little  con- 
fidence that  the  slope  is  different  from  zero  (or  other  values 
near  zero  that  might  be  thought  appropriate). 

How  much  of  the  scatter  in  the  source  estimates  is  real, 
that  is,  associated  with  real  differences  in  the  source  ampli- 
tude per  kilo ton  of  yield?  The  standard  deviations  in  Table 

4 are  about  35  to  45  percent  of  the  mean.  The  yield  is  only 
determined  to  within  about  t 15  percent  and  this  accounts 
for  some  of  the  scatter.  From  the  WWSSN  film  it  is  generally 
not  possible  to  measure  the  amplitudes  any  closer  than  within 

5 to  10  percent  (we  know  that  range  differences  add  errors 
of  a few  percent) . Also,  the  instrument  gain  may  vary  with 
time  and  not  be  the  same  as  the  assumed  values.  This  could 
add  errors  of  at  least  several  percent.  Therefore,  if  we 
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Figure  13.  The  4^  values  (at  0.02  KT)  are  plotted  versus  explosion 

yield  for  the  Yucca  Flat  and  Pahute  Mesa  events.  The  plot 
is  log-log  and  the  mean  and  95  percent  confidence  limits 
on  the  mean  are  shown.  One  decade  of  yield  is  shown  for 


TABLE  4 


INFERRED  VALUES  OF  f SCALED  TO  0.02  KT 


'f  (TUC) 


Events 

ALQ 

00 

TUC 

00 

(A: 

oo 

Yucca  Flat 

9 . 1 

13.2 

1.50 

95%  confidence  limits 

7.5 

- 10.9 

11.3  - 15.3 

1.  32 

Standard  deviation1 

44% 

35% 

24% 

n = 

17 

n = 18 

n = 

Pahute  Mesa 

6.6* 

11.2* 

1.57 

95%  confidence  limits 

4.8 

- 9.1 

7.8  - 16.0 

1.22 

Standard  deviation*" 

36% 

40% 

31% 

n = 

6 

n * 6 

n = 

PILEDRIVER 

3.9 

6.2 

1.59 

Deleting  one  event  known  to  be  above  the  water  table  and 
therefore  likely  to  have  lower  coupling. 


The  standard  deviation  is  given  as  a percentage  of  the 
mean. 
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assume  these  random  effects  are  independent  of  one  another, 
they  account  for  less  than  half  of  the  scatter.  We  also 
expect  some  real  source  amplitude  differences  due  to  variations 
in  ug  and  other  source  region  material  properties.  Any  non- 
axisymmetric  effects  that  are  present  may  also  cause  variations 
in  the  apparent  source  amplitude. 

In  summary,  we  have  assumed  the  source  was  a spherically 
symmetric  explosion  characterized  by  the  reduced  velocity 
potential  shown  in  Figure  9.  Correcting  for  the  burial  depth 
and  yield  according  to  (55,  we  estimated  the  required  to 
match  the  data  for  each  event.  Mean  values  were  computed  for 
the  Yucca  Flat  and  Pahute  Mesa  events  at  each  station.  The 
standard  deviations  are  about  35-45  percent  of  the  mean.  Scatter 
of  at  most  half  this  magnitude  is  expected  from  the  presence  of 
random  errors  in  the  yield  estimates  and  amplitude  measurements. 

Real  variations  in  the  for  each  population  are  therefore 
thought  to  have  standard  deviations  exceeding  20  percent 
of  the  mean  . 

There  are  several  features  of  the  '?  estimates  that  re- 

00 

quire  further  discussion.  First,  there  is  the  significance  of 

these  values  in  view  of  nonspherical  effects  such  as  tectonic 

release  and  spall  slapdown  that  are  widely  believed  to  influence 

Rayleigh  wave  generation.  Second,  there  is  the  asymmetry  in  the 

source  estimates  from  the  two  stations.  Finally,  there  are  the 

relative  values  for  the  three  test  areas.  The  ratio  of  the  mean 

Pahute  Mesa  'F  to  that  from  Yucca  Flat  is  0.73  at  ALQ  and  0.85 

at  TUC.  These  values  are  consistent, though  the  confidence 

limits  show  that  we  cannot  be  certain  that  the  lower  f for 

00 

Pahute  Mesa  is  meaningful.  Comparing  PILEDRIVER  to  Yucca  Flat, 
the  4^  ratios  are  0.43  for  ALQ  and  0.47  for  TUC.  Again,  the 
values  are  remarkably  consistent.  Within  the  limits  imposed 
by  the  statistics,  the  relative  source  levels  appear  to  be 
well  determined  by  the  procedure  we  have  used,  though  we  should 
keep  in  mind  that  it  is  the  relative  ug  'F08  that  is  being  determined. 
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The  absolute  source  levels  are  consistent  with  those 
computed  using  deterministic  models  for  the  constitutive  be- 
havior of  the  NTS  granite,  tuff  and  rhyolite.  For  example, 
in  Figure  14  we  show  computed  source  function  spectral  ampli- 
tudes, |$(a>)|,  for  several  NTS  materials.  These  source  func- 
tions are  from  Bache,  et  al. , (1975),  and  are  consistent  with 
the  relative  amplitudes  of  teleseismic  body  waves  which  are 
sensitiive  to  the  source  function  amplitude  near  1 Hz.  Varying 
the  air-filled  porosity  (<pQ)  of  the  Yucca  Flat  tuff,  we  get  the 
source  functions  shown  in  Figure  15.  This  figure  includes  the 
source  133  used  in  our  synthetic  seismogram  calculations.  The 
effect  on  the  Pahute  Mesa  rhyolite  source  of  varying  the  crush 
pressure,  P , is  shown  in  Figure  16.  The  P is  the  pressure 

G C 

at  which  the  air-filled  voids  are  irreversibly  removed  from 
the  rock. 

In  Figure  17  we  show  the  spectral  amplitude  of  the 
analytical  form  for  observed  source  functions  given  by  Haskell 
(1967).  The  "observed"  source  functions  were  computed  by 
Werth  and  Herbst  (1963)  from  ground  motion  measurements  in 
the  vicinity  of  several  low  yield  events.  Another  example 
is  given  in  Figure  18  where  the  PILEDRIVER  source  function 
130  is  compared  to  source  functions  computed  from  radial 
velocity  records  obtained  at  stations  at  two  horizontal  ranges 
(Perret,  1968) . 

The  theoretical  and  observed  values  from  Figures 
14  - 18  are  compared  to  those  inferred  in  this  study  in  Table 
6.  All  values  are  in  meters  cubed  and  scaled  to  0.02  KT. 
Clearly  our  inferred  values  are  about  the  right  size.  For 
Yucca  Flat  tuff  the  theoretical  source  with  low  (1.6  percent) 
air-filled  porosity,  source  133,  is  in  good  agreement  with 
the  observed  values.  The  questionable  value  for  the  observed 
tuff  V m is  from  the  RAINIER  event  and  the  material  properties 
were  probably  significantly  different  than  for  the  events  we 
studied.  For  Pahute  Mesa  the  theoretical  values  may  be  a 
little  small,  but  the  agreement  with  the  values  we  inferred 
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Frequency  (Hz) 


Figure  16.  Source  spectra  scaled  to  0.02  KT  are  shown  for 
calculations  124  (P  » 0.5)  and  125  (Pc  « 1.0) 
in  Pahute  Mesa  rhyolite. 


r 


Figure  17.  The  analytical  source  functions  calculated 
by  Haskell  (1967)  from  data  given  by  Werth 
and  Herbs t (1963)  are  shown  with  the 
amplitudes  scaled  to  0.02  KT  and  the  frequency 
axis  scaled  to  3 KT. 
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Ficjure  18.  The  source  function  amplitudes  computed  from  PILEDRIVER  observed 
velocity  records  at  two  raiujes  are  compared  to  the  theoretical 
source  130.  The  amplitude  axis  is  scaled  to  0.02  KT. 


is  quite  acceptable.  Finally,  for  PILEDRIVER  the  inferred 
values  are  significantly  smaller  than  the  observed  granite 
values  from  two  events  (Haskell's  value  is  for  HARDHAT)  or 
for  the  theoretical  value  which  was  computed  to  simultaneously 
match  the  near-field  and  teleseismic  short  period  data.  The 
reason  for  this  discrepancy  is  not  clear  at  this  time,  but  it 
is  not  obvious  which  values  are  more  nearly  correct. 

Another  interesting  comparison  is  to  the  values  in- 
ferred from  the  seismic  moments  computed  by  Tsai  and  Aki  (1971) 
and  summarized  by  Aki,  et  al.  (1974) . The  moment,  M , is  re- 
lated to  by 

Mq  = 47rpa2  . (6) 

The  P wave  velocity,  density  and  shear  modulus  for  our  deter- 
mination of  y are  (Figures  6 and  7)  : 

Yucca  Flat:  a = 2.35,  p = 1.86,  u = 31.4 

Pahute  Mesa:  a = 3.4  , p = 2.19,  yg  = 96.6 

PILEDRIVER:  a * 5.33,  p ■ 2.67,  y = 206. 

s 

Aki,  et  al.  (1974)  list  seismic  moments  for  eight  Yucca 
Flat  (MONERO,  BUFF,  PIRANHA,  BRONZE,  TAN,  OSCURO,  CORDUROY 
and  DUMONT)  and  four  Pahute  Mesa  (KNICKERBOCKER,  HALFBEAK, 
BENHAM  and  BOXCAR)  events  as  well  as  PILEDRIVER.  Using  (6) 
the  values  were  computed  for  each  of  these  events  and 
scaled  to  0.02  KT.  The  logarithmic  mean  values  for  each 
population  are  listed  in  Table  6 along  with  their  95%  con- 
fidence limits. 

The  moment  calculations  of  Tsai  and  Aki  (1971)  were 
based  on  ALQ  recordings  for  all  events  except  BOXCAR  and 
BENHAM  where  recordings  from  Spring  Hill,  Alabama  were  used. 

The  values  obtained  from  the  moment  are  slightly  smaller 
than  ours  for  PILEDRIVER  end  somewhat  larger  than  our  average 
for  the  other  two  populations.  In  view  of  the  differences 
in  the  methods  for  correcting  for  the  travel  path,  we  would 
not  expect  much  better  agreement. 
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Another  way  to  estimate  4^  was  suggested  by  Murphy 
(1974)  in  a discussion  of  the  paper  by  Aki,  £t  al.  (1974) . 

He  points  out  that  for  incompressible  materials 

V.  = R c/3  , (7) 

where  Rc  is  the  cavity  radius.  Using  an  empirical  relation- 
ship between  cavity  radius,  depth  and  yield  derived  by 
Orphal  (1970)  for  explosions  in  tuff,  he  finds  that 

= 1190  W0,89m3  . (8) 

The  values  computed  from  (8)  are  summarized  in  Table  6.  From 
a theoretical  point  of  view.  Equation  (7)  provides  an  upper 
limit  for  the  actual  'Vai.  For  the  calculations  leading  to  the 
theoretical  source  functions  shown  in  Figures  14  - 16,  the 

3 

ratio  of  R /3  to  V is  1.2  - 1.8  for  the  granite  and  rhyolites, 

5 to  14  for  the  saturated  tuffs  and  about  18  for  the  dry, 
porous  tuff.  Taking  these  factors  into  account,  the  values 
inferred  from  the  cavity  radii  suggest  that  our  inferred 
are,  perhaps,  a little  too  large. 

In  summary,  we  have  estimated  the  V x for  three  classes 
of  events,  assuming  the  source  is  spherically  symmetric.  For 
Yucca  Flat  and  Pahute  Mesa  the  theoretical  results,  which  were 
developed  to  match  the  teleseismic  short  period  body  waves, 
tend  to  support  the  smaller  of  the  two  estimates,  that  from 
the  ALQ  data.  The  values  inferred  from  the  cavity  radius  also 
support  the  ALQ  estimates.  (These  two  estimates  are  not 
entirely  independent  since  the  theoretical  calculations  took 
cavity  radius  into  account) . The  4^  values  inferred  from  the 
moment  are  larger  than  our  values,  but  this  is  likely  to  be 
due  to  the  different  path  correction  used  by  Tsai  and  Aki  (1971) . 
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TABLE  6 

COMPARISON  OF  THEORETICAL  AND  OBSERVED  ^ VALUES  TO  THOSE 
INFERRED  FROM  THE  ALQ  AND  TUC  RAYLEIGH  WAVES.  ALL 
VALUES  ARE  SCALED  TO  0.02  KT. 


ALQ  Inferred 


TUC  Inferred 


Inferred  from 
moment  values 
of  Aki,  et  al.  , 
1974 

Inferred  from 
cavity  radius 
(Murphy,  1974) 

Theoretical 

Figure  14 

Figures  15  & 16 

Observed 


Yucca  Flat 

9.1 

7.5  - 10.9 
n - 17 

13.2 

11.3  - 15.3 
n « 18 

16.6 

13.2  - 20.7 
n = 8 


15.1 

13.8  - 16.6 
n = 8 


4.8 

6.8,  10.3 


Figure  17  20.5* 

Figure  18 


Pahute  Mesa  PILEDRIVER 

6.9  3.9 

4.8  - 9.1 
n = 6 

11.2  6.2 

7.8  - 16.0 
n = 6 

16.2  3.75+ 

12.5  - 21.0 
n = 4 


12.2 

9.6  - 15.3 
n = 4 


4.2  9.2 

7.2 


10.0 

5.8,  16.7 


* 

Haskell  (1967)  lists  this  value  as  questionable. 
Based  on  a yield  of  56  KT. 
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For  PILEDRIVER  the  inferred  value  is  substantially 

smaller  than  the  theoretical  value  and  that  obtained  from 

integrated  velocity  recordings.  However,  in  this  case,  the 

estimate  from  the  moment  is  smaller  than  our  ¥ derived  from 

00 

either  ALQ  or  TUC. 

We  know  that  explosions  deviate  substantially  from 
spherical  symmetry.  In  the  next  section  we  study  the  in- 
fluence on  our  solution  of  some  plausible  nonspherical 
effects.  Another  important  aspect  of  our  solution  is  the 
difference  between  the  source  estimates  from  the  two  stations. 
This  will  be  discussed  in  Section  VII. 
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V.  CORRECTION  FOR  THE  EFFECT  OF  A 
DOUBLE -COUPLE  ENHANCEMENT  OF  THE  SOURCE 


Two  often  observed  and  discussed  phenomena  that  are 
incompatible  with  the  assumption  of  spherical  symmetry  for 
underground  nuclear  explosions  are  surface  spall  and  tectonic 
strain  release.  In  the  next  two  sections  we  estimate  the 
effect  these  phenomena  might  have  on  our  source  estimates. 

The  evidence  for  explosions  being  accompanied  by  a 
double-couple  generated  by  some  type  of  tectonic  strain  re- 
lease is  quite  convincing  (Toksoz,  et  al. , 1971;  Lambert,  et 
al. , 1972;  Aki  and  Tsai,  1972).  Toksoz  and  Kehrer  (1972) 
inferred  the  size  of  the  double-couple  associated  with  several 
of  the  events  studied  here.  Assuming  a vertical  strike-slip 
orientation,  the  Rayleigh  waves  from  the  composite  source 
made  up  of  the  explosion  plus  the  double-couple  may  be  written 

w * w (1  + F sin  29)  , (9) 

where  wfc  and  wg  are  the  total  and  explosion  generated  Rayleigh 
wave  amplitudes.  The  F represents  the  size  of  the  double- 
couple relative  to  the  explosion  and  9 is  the  azimuth  from  the 
strike.  Using  the  F and  9 given  by  Toksoz  and  Kehrer  (1972), 
the  wfc  were  computed  for  several  events  and  are  tabulated  in 
Table  7.  Also  tabulated  are  the  TUC/ALQ  ratios.  These  in- 
dicate the  relative  size  of  the  source  as  viewed  from  TUC  com- 
pared to  its  size  when  viewed  from  ALQ.  For  the  Yucca  Flat 
and  Pahute  Mesa  events  the  double-couple  radiation  pattern 
tends  to  make  the  source  look  bigger  at  TUC.  If  the  population 
of  events  in  Table  7 is  a random  sample,  we  could  conclude 
that  the  source  appears  to  be  something  like  5 to  30  percent 
larger  at  TUC  due  to  the  tectonic  release  radiation  pattern. 
This  would  explain  part  of  the  difference  in  the  fa5  estimates 
from  ALQ  and  TUC  in  Table  4.  However,  the  inferred  radiation 
pattern  for  PILEDRIVER  tends  to  slightly  increase  the  ALQ-TUC 
asymmetry. 
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TABLE  7 


CORRECTION  OF  THE  RAYLEIGH  WAVE  AMPLITUDES  FOR  THE  DOUBLE -COUPLE 


SOLUTION  OF  TOKSOZ 


EVENT 

F 

0 

Yucca  Flat 

CUP 

.55 

ro 

o 

0 

TAN 

.39 

347° 

BRONZE 

.33 

5° 

BUFF 

.31 

28° 

CORDUROY 

.72 

347° 

BILBY* 

.47 

340° 

HAYMAKER* 

.33 

340° 

Mean  + one  standard 

deviation 

Pahute  Mesa 

DURYEA 

.75 

355° 

CHARTREUSE* 

.90 

353° 

HALFBEAK* 

.67 

345° 

BENHAM* 

.85 

345° 

BOXCAR* 

.59 

346° 

GREELEY* 

1.6 

355’ 

Mean  + one  standard  deviation 


PILEDRIVER  3.2  340° 


AND  KEHRER  (1972) 


TUC 

wt/we 

ALQ 

V”e 

(TUC) 

wfc(ALQ) 

1.43 

.87 

1.64 

1.34 

1.31 

1.02 

1.33 

1.09 

1.22 

1.18 

.76 

1.55 

1.64 

1.57 

1.05 

1.35 

1.43 

.94 

1.25 

1.30 

.96 

36+ . 15 

1.19+.30 

1.20+.29 

1.74 

1.44 

1.21 

1.87 

1.58 

1.18 

1.58 

1.56 

1.01 

1.74 

1.  71 

1.02 

1.52 

1.48 

1.03 

2.58 

1.94 

1.33 

84+ . 38 

1 . 62+ . 18 

1.13+.13 

3.00+ 

3.55+ 

.85 

‘Events  not  included  in  the  population  observed  at  ALQ  and/or 
TUC 

+These  values  are  based  on  synthetic  seismogram  calculations. 
They  are  about  10  percent  smaller  than  those  obtained  from 
Equation  (9)  due  mainly  to  the  source  amplitude  being  slightly 
different  at  the  measurement  period  than  at  zero  frequency. 
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From  the  values  in  Table  7 we  conclude  that  the  double- 
couple component  increases  the  Rayleigh  waves  at  ALQ  and  TUC 
in  most  cases.  If  we  correct  for  this  effect,  our  estimate 
of  the  ^ for  the  explosion  is  reduced  accordingly.  We  should 
also  point  out  that  for  a double-couple  of  this  size  and 
orientation  the  group  delay  at  the  source  is  quite  small.  As 
long  as  F < 1,  the  double-couple  has  almost  no  affect  on  the 
seismogram  except  to  scale  its  amplitude.  The  same  is  true 
for  larger  F if  we  are  on  a positive  lobe  of  the  radiation 
pattern.  This  is  illustrated  in  Figure  19  where  the  theoretical 
seismogram  for  PILEDRIVER  at  TUC  from  Figure  5 is  compared  to 
similar  seismograms  with  the  source  including  a strike-slip 
double-couple  component. 


The  double-couple  components  for  the  seismograms  of 
Figure  19  are  computed  using  the  tectonic  release  model  of 
Archambeau  (1972).  If  the  tectonic  stress  release  has  the 
sense  of  pure  strike-slip  faulting,  the  F factor  is  related 
to  the  parameters  of  the  model  by 


F = if 


qxy  Ro 


where  axy  is  the  horizontal  stress  drop  and  Rq  is  the  radius 

of  the  spherical  volume  in  which  this  (average)  stress  drop 

occurs.  For  example,  if  the  PILEDRIVER  yield  is  56  KT  and 

u * 206  kbar,  several  values  of  R , a and  V at  0.02  KT  • 
s o xy  00 

that  lead  to  F * 3.2  are  given  below 


f«02  (m3)  Ro  (m)  qxv  (bars) 


4 

609 

100 

4 

484 

200 

2 

300 

419 

8 

609 

200 

0 
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Interpretation  of  F in  terms  of  a plane  strike-slip 
dislocation  can  also  be  made.  The  moment  for  the  Archambeau 
source  is  (Minster  and  Suteau,  1977) 


°xy  Ro  ' 


while  for  a dislocation  source  it  is 

M = u S D 
o s 

where  S is  the  fault  area  and  D is  the  average  dislocation, 
Therefore, 


F = 


S D 

a-rr  v 


and  a table  like  that  given  above  for  the  Archambeau  source 
is  easily  constructed. 

In  Table  8 are  given  the  corrected  values  of  com- 
puted by  dividing  the  mean  values  of  Table  4 by  the  mean 
double-couple  enhancement  from  Table  7.  These  corrections 
are  rather  crude,  but  do  indicate  the  trend  of  the  correction 
and,  roughly,  its  magnitude.  For  Yucca  Flat  and  Pahute  Mesa 
the  corrected  solutions  are  quite  satisfactory  and  might  be 
preferred  to  those  in  Table  4.  For  PILEDRIVER  the  new 
values  in  Table  8 are  quite  small  and  are  difficult  to  justify 
when  compared  to  theoretical  and  empirical  data  (see  Table  6) . 
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50  sec 


Figure  19.  Theoretical  seismograms  for  PILEDRIVER  at  TUC  are 

shown  for  a composite  source  including  the  explosion 
plus  a strike-slip  double-couple.  The  azimuth 
is  157®  clockwise  from  the  strike  and  the  relative 
size  of  the  double-couple  is  indicated  by  the  F 
factor  with  a negative  F signifying  that  the  station 
is  on  the  negative  lobe. 
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TABLE  8 

INFERRED  V VALUES  CORRECTED  BY  THE  DOUBLE -COUPLE 
00 

FACTORS  FROM  TABLE  7 


y«(Tuc) 


Events 

ALQ 

TUC  ¥ 

00 

CALQ) 

Yucca  Flat 

7.6 

9.7 

1.28 

Pahute  Mesa 

4.1 

6.1 

1.49 

PILEDRIVER 

1.1 

2.1 

1.88 
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VI.  CORRECTION  FOR  THE  EFFECT  OF 
SURFACE  SPALLATION 


6 . 1 Introduction 

Underground  nuclear  explosions  almost  always  cause 
spallation  near  the  free  surface  (e.g.,  Eisler  and  Chilton, 

1964)  . When  the  spalled  material  falls  back  to  the  surface 
(spall  closure) , some  impulse  is  delivered  and  a large  signal 
associated  with  spall  closure  can  often  be  observed  on  near- 
field recordings  (e.g.,  Eisler,  et  al. , 1966;  Toman,  et  al. , 

1973) . These  data  can  be  used  to  estimate  the  dynamics  of  the 
spall  zone  and  its  extent  (e.g.,  Chilton,  et  al. , 1966;  Viecelli, 
1973;  Sobel,  1978).  Viecelli  (1973)  suggested  that  the  spall 
closure  makes  a large  contribution  to  the  recorded  Rayleigh 
waves  in  the  far-field.  If  this  were  true,  it  would  strongly 
compromise  the  results  of  our  analysis  in  previous  sections. 

In  this  section  we  first  summarize  the  estimates  for  the 
extent  of  spall  given  by  Viecelli  and  by  Sobel  (1978 ) . Following 
Viecelli,  we  use  these  estimates  to  bound  the  amplitude  of  the 
impulse  applied  by  spall  closure.  Using  some  theoretical 
results  given  by  Harkrider,  et  al. , (1974),  we  then  compute 
the  Rayleigh  waves  due  to  spall. 

If  large  amounts  of  energy  are  trapped  in  the  spall 
zone,  we  should  also  account  for  the  loss  of  this  energy  from 
the  calculation  of  the  explosion  produced  Rayleigh  waves.  We 
approximate  this  by  deleting  from  the  Rayleigh  wave  calculation 
all  or  a portion  of  the  upward  traveling  energy  from  the  ex- 
plosion. The  Rayleigh  waves  generated  by  the  spall  closure 
are  then  superimposed  on  this  truncated  source. 

In  summary,  to  account  for  the  effect  of  surface 
spallation  it  is  necessary  to  subtract  some  portion  of  the 
explosion  energy  and  then  to  add  the  energy  from  the  spall 
closure.  In  view  of  the  approximations  involved,  this  cannot 
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be  done  with  much  certainty.  However,  experiments  with 
plausible  values  for  the  governing  parameters  allow  estimation 
of  some  bounds  on  the  spall  contribution. 


I 
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6 . 2 Estimation  of  the  Impulse  due  to  Spall  Closure 

Using  data  from  surface  accelerometers,  some  estimates 
can  be  made  for  the  extent  of  the  spall  region  (Eisler  and 
Chilton,  1964;  Viecelli,  1973;  Sobel,  1978).  Viecelli  (1973) 
studied  the  data  for  six  explosions,  two  in  Rainier  Mesa  tuff 
and  four  in  Pahute  Mesa  tuff  and  rhyolite.  Assuming  a density 
of  2 gm/cm2,  he  estimated  the  mass  of  the  upthrown  material 
to  be 


<M>  s*?  1.6  x 1012  W grams  , (10) 

where  W is  the  yield  in  kilotons.  Sobel  (1978)  carried  out 
a similar  study,  looking  at  the  near-field  data  from  some 
twenty-six  events  in  a variety  of  materials.  Using  her  derived 
relations  for  the  dimensions  of  the  spall  zone  and  p = 2 gm/cm2, 
the  estimate  for  the  mass  is 

<M>  =s  10  x 1012  W97  grams  . (11) 


The  factor  of  six  difference  between  the  two  estimates  is 
mainly  due  to  different  methods  for  estimating  the  maximum 
extent  of  the  spall  zone. 

Assuming  a conical  displacement  profile,  Viecelli 
estimated  the  total  downward  impulse  due  to  the  spall  plate 
falling  back  to  be 


<I> 


1/2 

( 2g  V 


(12) 


y 
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where  g is  the  gravitational  acceleration  and  hQ  is  the  maximum 
height  of  spall.  Using  Viecelli's  estimates  for  the  spall 
dimensions, 

<I>  as  4.6  x 10 ^4  W dyne-sec;  (13) 

while  the  dimensions  found  by  Sobel  lead  to 

<I>  as  8.6  x 1014  W1*15  dyne-sec  . (14) 

At  100  KT  these  estimates  differ  by  a factor  of  3.7. 

The  temporal  and  spatial  distribution  of  the  spall 
closure  is  quite  small  compared  to  the  periods  and  wave- 
lengths of  primary  interest  and  it  seems  reasonable  to  represent 
this  source  as  an  impulse  load  at  the  surface.  The  question 
is  then,  how  good  are  the  estimates  in  (13)  and  (14)  for  the 
amplitude  of  this  impulse? 

The  mass  was  computed  by  assuming  that  the  volume  of 

the  spall  can  be  approximated  by  a circular  disk  of  uniform 

2 1/2 

thickness.  Then  the  impulse  is  proportional  to  org  D hQ' 
where  p is  density,  D is  the  thickness  of  the  spall  plate 

and  r is  its  radius.  The  spall  thickness  is  about  100m 

s 1/3 

for  100KT  and  is  proportional  to  W ' . At  depths  this  shallow 

the  average  density  is  probably  closer  to  1.7  gm/cm^  than  to 

2 gm/cm^ , the  value  used  to  compute  (10)  and  (11).  Using 

- p = 1.7,  our  impulse  estimates  are  reduced  by  15  percent. 

The  hQ  is  determined  from  surface  measurements  and 
is  attributed  entirely  to  separation  at  the  single  depth  D. 

Since  the  total  surface  displacement  includes  contributions 
from  spallation  at  many  depths  (Eisler,  et  al. , 1966)  the 
latter  assumption  is  likely  to  lead  to  an  overestimate  of  the 
spall  impulse,  even  if  the  estimates  for  rg  and  D are  accurate. 


■BMggasa  aaaapfcaa&aac 
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If  we  assume  N spall  plates  of  thickness  D ^ separated  by' 
openings  hj,  j = 1,...,N,  then  the  total  impulse  is  propor- 
tional to 


t 


1 E hi 


(15) 


\l/2 

Ji  V ^ 7 ‘ 

N N 

Given  that  D.  = D and  . h.  = h , it  is  easily  shown  that 

1=1  1 1=1  1 ° 1/2 


this  factor  must  be  less  than  or  equal  to  Dh 


its  N— 1 value. 


The  conclusion  that  the  impulse  in  (14)  is  an  over- 
estimate was  also  drawn  by  Sobel  (1978)  after  studying  sub- 
surface accelerometer  records.  She  found  the  spall  energy 
leaving  the  source  region  to  be  much  smaller  than  expected 
from  her  estimates  of  the  extent  of  spall  and  attributed  the 
discrepancy  to  the  pressure  of  multiple  spall  openings. 

Another  question  to  be  addressed  is  the  time  lag 
between  the  explosion  and  the  spall  closure.  It  is  easily 
shown  that  this  time  lag  TL  is 


TL  = Ts  + 2 v/2g  hQ. 


(16) 


where  T is  the  time  when  spall  parting  occurs.  For  explosions 


,1/3 


1/3 


in  tuff  at  a scaled  depth  of  122  W , Tg  % .06  W . At 
100  KT  this  is  0.28  seconds.  The  second  term  in  (16)  is  the 
time  the  spall  remains  open.  If  there  is  only  one  spall 
opening,  Sobel 's  results  give  .52  for  this  time  (1.2 

seconds  at  100  KT)  while  Viecelli  finds  a yield  independent 
value  of  1.75  seconds.  Therefore,  for  yields  of  50  to  200 
KT  the  total  T-  is  1.3  to  2.0  seconds,  assuming  a single 

Li 

spall  parting.  For  multiple  spall  openings  the  second  term 
in  (16)  is  reduced  the  same  way  the  impulse  is  reduced  by 
multiple  spall  partings.  Since  this  term  dominates  and  multiple 
spall  partings  do  occur,  a lag  of  about  1 second  is  thought 
to  be  most  appropriate. 
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We  assume  that  spall  closure  can  be  represented  as  an 
applied  normal  force  at  the  surface  of  a multilayered  halfspace. 
Harkrider,  et  al. , (1974)  gives  the  expressions  for  the  Rayleigh 
waves  due  to  loading  of  this  kind.  The  vertical  (positive 
downward)  displacement  is 

ws(r,o))  = -iH  Ar  P (oj)  ^t)  ' d7) 

where  we  have  introduced  the  notation 


A„  = -i  <GN  - LH)  (18) 

for  the  amplitude  excitation  following  Harkrider  (1964).  The 
source  term  is 

CO 

P(u)  = j Po(r,u)  JQ{~,  r dr  , (19) 

where  P (r,u))  is  an  azimuthallv  symmetric  Fourier  time- 
o 

transformed  normal  stress  distributed  on  the  surface.  If  the 
spatial  extent  and  time  duration  of  the  surface  load  are 
small  compared  to  wavelengths  and  periods  of  interest,  the 
total  impulse  is 

<I>  = 2H  P . (20) 

Using  separate  models  for  the  local  source  region 
and  the  remainder  of  the  path,  the  expression  for  the  spall 
induced  Rayleigh  waves  that  is  analogous  to  Equation  (1)  tor 
an  explosion  is 


ws(r,uO  = -in  P(cij)  Ar 


T(uj)  HC2)(  £!£  ) e'yr  . 
o 


l r \ 

y ae  sm  A j 


(21) 


t( 


■1 


I) 


I 


where  terms  for  the  attenuation  and  sphericity  correction  have 
also  been  introduced. 

Comparing  (1)  and  (21)  we  see  that  a spherically 
symmetric  explosion  and  a spall  closure  source  give  identical 
Rayleigh  waves  when 

/S  A 

4Ms  » ^ = i P . (22) 


For  explosions  in  Yucca  Flat  tuff  our  models  give  'Ks  /c^j  to 
be  about  0.39  sec/km  at  10  seconds.  Then  the  spectral  amplitudes 
"for  the  two  sources  are  the  same  when 

ll  02 

2.45  x 10XJ  W = P , (23) 

02  ? 3 

where  is  the  zero  frequency  value  of  1 in  m scaled  to 

C.02  KT.  If  the  scaled  V x is  about  10,  a reasonable  value 

for  Yucca  Flat  tuff  (Table  4) , the  impulse  calculated  by 

Viecelli,  Equation  (13),  gives  Rayleigh  waves  about  30  percent 

• • 

as  large  as  those  from  the  explosion.  For  yields  around  100  KT 
Sobel's  estimate,  Equation  (14),  gives  Rayleigh  waves  that 
are  the  same  size  as  those  from  the  explosion.  We  have  already 
said  that  the  impulse  values  in  Equations  (13)  and  (14)  appear 
to  be  overestimated.  However,  they  are  large  enough  that  we 
must  explore  their  potential  influence  more  carefully. 
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In  Figure  20  we  show  the  effect  of  adding  Rayleigh 
waves  from  the  spall  closure  impulse  to  the  Rayleigh  waves 
from  the  spherically  symmetric  source  for  a typical  Yucca  Flat 
event.  The  top  two  records  are  the  TAN  observations  and 
synthetic  seismograms  from  Figure  5.  We  then  show  the  Rayleigh 
waves  from  the  spall  closure  computed  using  Equation  (21). 

The  PU)  was  assumed  to  be  an  impulse  source.  Equation  (20), 
so  the  amplitude  of  the  signal  is  directly  proportional  to  the 
size  of  the  impulse.  The  last  three  records  show  the  effect 
of  summing  the  explosion  and  spall  produced  Rayleigh  waves. 
There  are  two  parameters  for  the  summation,  the  lag  time  Tt 
and  the  ratio  of  impulse  and  explosion  sources,  P /Y  , which 
is  motivated  by  (23)  . In  this  context  '?  = 4''02  W/0.02. 

The  summed  seismograms  in  Figure  20  were  done  for  two 
lag  times  and  two  P/'C  ratios.  For  an  event  like  TAN,  a one 
second  Tt  seems  more  appropriate,  though  somewhat  larger  values 
are  possible.  The  P/C^  ratios  can  be  interpreted  in  terms 
of  the  impulse  and  compared  to  the  estimates  of  Viecelli  and 
Sobel  in  Equations  (13)  and  (14) . From  (20)  we  have 

<I>  - 211  n ( 6 . 2 x 1012)  t;02  W , (24) 

where  n is  1 or  2 for  the  seismograms  in  Figure  20.  Then  if 
4^  2 = 11.8  m , the  impulse  is  1.0  or  2.0  times  the  impulse 
estimate  of  Viecelli.  For  n * 2,  the  impulse  is  1.7  to  2.1  . 

times  smaller  than  Sobel* s estimate  at  50  to  200  KT.  It  is 
probably  not  reasonable  to  expect  an  impulse  much  larger  than 
this. 

The  first  thing  we  notice  about  the  composite  records 
in  Figure  20  is  that  the  waveform  is  not  greatly  changed  by 
the  addition  of  the  spall  closure  contribution.  In  fact,  at 
ALQ  the  change  is  hardly  perceptible.  At  TUC  the  second 
trough  is  reduced  and  the  agreement  with  the  observed  record 
is  actually  enhanced.  The  relative  amplitudes  of  the  phase 
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Figure  20.  The  effect  of  the  spall  closure  impulse  on  the  Rayleigh  waves  for 
Yucca  Flat  explosions  is  shown  for  two  lag  times  and  two  spall 
impulse  amplitudes. 


corresponding  to  the  maximum  peak-to-peak  excursion  on  the 
observation  is  also  indicated  on  the  plots.  The  effect  is 
a little  greater  at  TUC  than  at  ALQ,  but  is  less  than  20 
percent  for  every  case. 

Since  the  effect  of  the  spall  impulse  contribution 

is  to  reduce  the  Rayleigh  wave  amplitude,  we  would  account 

for  this  effect  by  increasing  our  estimate  for  the  'i  due 

00 

to  the  explosion  itself_.  _The  increase  is,  of  course,  by  the 
inverse  of  the  relative  amplitudes.  The  new  inferred  values 
are  shown  on  Figure  20.  (The  spall  closure  impulse  is  also 
increased  by  the  same  amount  since  P/'F^  is  fixed)  . These 
values  are  not  specifically  for  TAN  but  are  based  on  the  mean 
Yucca  Flat  'Foo  given  in  Table  4.  The  effect  on  the  TUC/ALQ 
ratio  is  also  shown  on  the  plot.  The  addition  of  a spall  con- 
tribution causes  the  discrepancy  between  the  ALQ  and  TUC  esti- 
mates to  be  reduced  a little. 

The  spall  contribution  might  be  expected  to  be  quite 
variable  and  could  account  for  some  of  the  scatter  in  our 
estimates  for  Yucca  Flat  events,  though  it  would  be  sur- 
prising if  it  introduced  scatter  much  bigger  than  10  to  20 
percent.  Some  evidence  for  the  consistency  of  this  contribution 
is  the  remarkable  similarity  in  the  waveforms  for  different 
events  (Figure  1)  . 

What  about  the  effect  of  spall  closure  on  the  signals 
from  events  in  other  materials,  the  Pahute  Mesa  tuffs  and 
rhyolites  and  the  PILEDRIVER  granite?  We  have  no  reason  to 
suppose  the  dimensions  of  the  spall  zone  are  dependent  on  the 
source  material.  Viecelli's  estimates  are  actually  based 
predominantly  on  observations  of  Pahute  Mesa  events.  Sobel 
looked  at  events  in  many  materials  but  no  systematic  dependence 
on  source  material  is  discernible.  Since  the  density  of  the 
near  surface  materials  is  likely  to  be  pretty  much  the  same 
at  all  these  areas,  we  can  only  conclude  that  the  impulse 
estimates  in  (13)  and  (14)  are  valid  for  all  three  areas. 
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The  major  factor  that  distinguishes  Yucca  Flat  from  the 
other  testing  areas  in  regard  to  the  contribution  from  spall 
closure  is  the  substantially  smaller  value  of  u 4^  for  the 
Yucca  Flat  events.  In  view  of  (22),  this  means  that  a larger 
impulse  is  required  to  have  the  ratio  of  the  spall  closure 
to  explosion  Rayleigh  waves  be  the  same  at  Pahute  Mesa  and 
PILEDRIVER  as  it  is  at  Yucca  Flat.  The  pertinent  numbers  are 
given  in  Table  9 . _ We  see  that  the  impulse  must  be  about  twice 
as  big  for  Pahute  Mesa  and  PILEDRIVER  events  to  have  the  same 
relative  contribution  as  at  Yucca  Flat.  From  the  Yucca  Fiat 
example  in  Figure  20,  one  would  conclude  that  the  potential 
contribution  of  the  spall  closure  can  be  neglected  for  the 
other  two  areas.  However,  before  closing  this  subject,  we 
need  to  look  at  the  effect  of  deleting  the  portion  of  the 
explosion  energy  that  is  absorbed  in  the  spallation  process. 


6 . 4 Modification  of  the  Explosion  Source  to  Account  for 
Energy  Lost  to  Spallation 

It  is  inconsistent  to  add  a downward  impulse  source 
to  represent  the  spall  closure  while  continuing  to  model  the 
explosion  as  a spherically  symmetric  point  source  in  a layered 
elastic  half space.  The  energy  for  the  spallation  comes  from 
the  waves  propagated  upward  from  the  source.  Thus,  if  a large 
amount  of  spallation  occurs,  we  expect  a substantial  reduction 
in  the  body  wave  phase  pP  from  what  is  predicted  by  elastic 
theory . 

Our  theoretical  formulation  for  computing  Rayleigh 

waves  can  be  modified  to  allow  total  or  partial  suppression 

of  the  source  radiation  emitted  upward  or  downward  from  the 

source.  For  a spherically  symmetric  source  this  is  done  by 

replacing  the  depth-dependent  excitation  factor  K in 

S1 

Equation  (11)  by 
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r = 1:  All  upgoing  waves  suppressed, 

F * 0:  Total  solution 

T = -1:  All  downgoing  waves  suppressed. 

The  stress  and  displacement  eigenfunctions  in  the  definitions 
of  Ks^,  Ms ^ are  defined  by  Harkrider  (1964,1970)  . The  T can 
take  any  value  in  the  range  [-1,1],  and  represents  the  propor- 
tion of  up  or  downgoing  waves  suppressed. 

In  the  formulation  of  Equation  (25)  F may  be 
constant  or  it  may  be  frequency  dependent.  In  the  applications 
here,  we  will  assume  a constant  F,  mainly  because  we  have  no 
firm  basis  for  doing  otherwise.  As  far  as  deleting  the  energy 
lost  to  spallation  is  concerned,  it  might  be  more  reasonable 
to  suppose  that  this  occurs  mostly  at  the  high  frequencies 
and  the  low  frequencies  are  much  less  affected.  More  theoret- 
ical work  is  required  to  develop  this  line  of  reasoning  much 
farther. 
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In  a previous  report  (Bache,  Goupillaud  and  Mason, 

1977) , we  studied  the  effect  of  suppressing  various  proportions 
of  the  upgoing  waves  on  the  Rayleigh  waves  at  3,000  km.  The 
two  cases  examined  were  like  the  Yucca  Flat  and  PILEDRIVER 
events  studied  here.  For  the  Yucca  Flat-like  event  we  found 
that  suppression  of  the  upgoing  wave  tends  to  increase  M 
compared  to  the  total  source  because  the  upgoing  and  downgoing 
portions  of  the  wave  field  generate  Rayleigh  waves  that 
destructively  interfere.  On  the  other  hand,  for  the  PILEDRIVER- 
like  event,  suppression  of  the  upgoing  waves  substantially 
reduced  Mg  because  the  Rayleigh  waves  from  the  two  parts  of 
the  source  are  in  phase  and  that  from  the  upgoing  waves  is 
considerably  larger.  These  conclusions  were  based  on  the 
well  dispersed  20  second  waves  on  the  teleseismic  records. 

However,  similar  results  are  obtained  for  the  Airy  phase  of 
interest  at  ALQ  and  TUC. 

The  different  contributions  of  the  upgoing  and  downgoing 
waves  for  events  in  granite  compared  to  those  in  Yucca  Flat 
tuff  can  be,  in  large  part,  attributed  to  the  importance  of 
S waves  for  Rayleigh  wave  generation.  For  spherically  symmetric 
explosions  in  granite  the  only  large  source  of  S waves  is  the 
pS  conversion  at  the  free  surface  which  is  deleted  when  the 
upgoing  waves  are  suppressed.  However,  there  is  a strong 
impedence  contrast  below  the  source  for  Yucca  Flat  events 
(Figure  6) , which  continues  to  generate  S waves  in  the  absence 
of  the  pS  conversion  at  the  free  surface. 

In  Figure  21  we  show  the  following  seismograms  for 
our  representative  Yucca  Flat  event  (TAN)  at  ALQ  and  TUC:  j 

A.  The  observation.  | 

B.  The  synthetic  for  the  spherically  symmetric  \ 

explosion  (Figures  5 and  20) . I 

t 

i 

C.  The  synthetic  for  the  source  with  all  upgoing  i 

waves  suppressed.  j 
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D.  A composite  synthetic  including  the  total  source 

plus  the  spall  contribution  with  = 1.24  x 101 

and  Tl  = 1 second  (Figure  20) . 

E.  A composite  synthetic  with  only  the  downgoing  waves 
from  the  explosion  (C)  and  the  same  spall  contri- 
bution as  in  D. 

F . A composite  synthetic  with  half  the  upgoing  waves 
suppressed  and  the  same  spall  impulse  as  D. 

G.  A composite  synthetic  with  25  percent  of  the 
upgoing  waves  suppressed  and  the  same  spall 
impulse  as  D. 

A key  feature  of  these  synthetic  seismograms  is  that 
the  amplitudes  vary  within  about  ±10  percent  of  that  for  the 
spherically  symmetric  point  source  model.  Thus,  our  assumptions 
about  the  spalling  and  degradation  of  the  upgoing  waves  are 
not  too  important.  For  each  seismogram  in  Figure  21,  we  also 
list  the  inverse  of  the  relative  amplitude  times  the  mean 
Yucca  Flat  fm  from  Table  4.  This  shows  the  amount  by  which 
the  inferred  values  are  corrected  by  the  effects  represented 
in  each  synthetic  record. 

The  span  of  reasonable  models  is  believed  to  be 
represented  by  seismograms  B,  C,  F and  G,  though  the  impulse 
contribution  may  be  overestimated.  (For  = 11.8,  the 
impulse  has  the  amplitude  estimated  by  Viecelli,  Equation 
(13)).  These  indicate  that  the  effect  of  spall-related 
phenomena  is  likely  to  increase  the  inferred  by  a small 
amount,  perhaps  5 to  10  percent.  The  discrepancy  between 
the  ALQ  and  TUC  estimates  is  also  reduced  slightly. 

Some  evidence  for  the  preferred  model  can  be  deduced 
from  the  comparison  of  theoretical  and  observed  waveforms. 

The  explosion  itself,  Case  A,  provides  a good  fit  at  both 
stations  and  addition  of  the  spall  closure  impulse,  Case  D, 
seems  to  improve  it.  However,  the  fit  is  substantially 
degraded  by  deleting  the  upgoing  waves.  Therefore,  the  best 
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fit  to  the  waveform  is  by  the  seismograms  in  which  little 
or  none  of  the  upgoing  waves  are  deleted,  Cases  D and  G. 
Taking  all  factors  into  consideration,  we  prefer  the  model 
for  Case  G. 


For  PILEDRIVER,  or  other  explosions  in  granite,  the 
story  is  quite  different.  Our  assumptions  about  spallation, 
and  especially  the  degradation  of  the  upgoing  waves,  have 
a strong  influence  on  the  amplitude  of  the  synthetic  seismo- 
grams. Synthetic  seismograms  for  granite  are  shown  in 
Figure  22.  The  seismograms  are  analogous  to  those  in  Figure 
21  except  that  the  amplitude  of  the  spall  closure  impulse 
has  been  adjusted  to  be  about  the  same  for  Cases  D to  G while 
the  inferred  varies  over  a wide  range.  The  impulse  values 
for  the  mean  of  the  ALQ  and  TUC  inferred  are  as  follows: 


Case 

' rri  ' 

Inferred 

^ • 0 2 
* 00 

<I> 

( dyne- 

sec) 

D 

4.1  x 1011 

5.3 

6.8 

X 

1014 

w 

E 

1.24  x 1011 

19.6 

7.6 

X 

io14 

w 

F 

2.1  x 1011 

8.3 

5.3 

X 

1014 

w 

G 

2.1  x 1011 

6.4 

4.1 

X 

1014 

w 

Therefore,  the  impulse  added  to  the  source  is  about  the  same 
size  as  for  the  Yucca  Flat  event  in  Figure  21  and  is  comparable 
to  the  value  (4.6  x 10^4  W)  found  by  Viecelli. 

The  major  factor  controlling  the  PILEDRIVER  seismograms 
is  the  proportion  of  upgoing  waves  deleted.  The  amplitude 
of  the  seismogram  for  the  total  source  is  nearly  four  times 
as  large  as  that  from  the  downgoing  waves  alone  (compare  B and 
C) . The  spall  closure  impulse,  as  expected,  has  a relatively 
small  effect.  Therefore,  the  amount  of  upgoing  energy  that 
is  lost  to  spall  or  scattered  is  of  great  importance  for  deter- 
mining the  granite  source.  The  waveform  fit  does  not  help 
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us  to  choose  between  the  alternatives.  The  fit  to  the  ALQ 
observations  is  equally  good  for  all  cases.  At  TUC  the  fit 
is  poorer  than  at  ALQ  and  is  not  much  better  for  one  case 
than  for  another. 

The  analysis  of  PILEDRIVER  is  further  complicated  by 
the  presence  of  a large  double-couple  component  in  the  source 
for  surface  waves.  As  we  pointed  out  in  Section  V,  the 
double-couple  solution  of  Toksoz  and  Kehrer  (1972)  predicts 
double-couple  radiation  that  dominates  the  Rayleigh  wave 
(Table  7) . The  waveform  from  the  double-couple  is  nearly 
the  same  as  that  from  the  explosion  alone  (Figure  19) . We 
must  take  all  these  factors  into  account  when  developing  our 
final  solution  for  PILEDRIVER. 

The  events  at  Pahute  Mesa  are  more  nearly  like  PILEDRIVER 
than  the  Yucca  Flat  events  as  far  as  the  effect  of  spallation 
is  concerned.  The  pertinent  synthetic  seismograms  are  shown 
in  Figure  23  in  the  same  format  as  used  in  Figure  21  and  22, 
except  that  Cases  C and  E are  not  shown.  The  amplitude  of  the 
spall  closure  impulse  for  these  seismograms  is  as  follows: 


Case 

p/T.(drne:sec) 

Inferred 

<I> 

(dyne-sec) 

D 

1.86  x 1011 

9.5 

5.6 

x 1014  W 

F 

2.48  x 1011 

11.1 

8.6 

x 1014  W 

G 

2.48  x 1011 

13.0 

10.1 

x 1014  W 

This  impulse  ranges  from  1.2  to  2.2  times  as  large  as  that 
found  by  Viecelli  (1973)  for  Pahute  Mesa  events. 

As  with  PILEDRIVER,  suppression  of  a portion  of  the 
upgoing  waves  has  a strong  effect  on  the  amplitude  of  the 
explosion  source  required  to  match  the  observations.  We  also 
note  that  the  spallation  effects  seem  to  increase  the  difference 
between  the  ALQ  and  TUC  estimates.  The  agreement  of  observed 
and  theoretical  seismograms  is  quite  good  for  both  stations 
and  seems  to  be  independent  of  our  assumptions  about  the  source. 
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Figure  23.  Synthetic  seismograms  for  a variety  of  source  representations 
for  a Pahute  Mesa  event. 


An  important  conclusion  of  this  analysis  of  the  spall 
effects  is  that  it  is  not  the  spall  closure  impulse  that  is 
important,  but  the  possible  loss  of  energy  from  the  upgoing 
waves.  At  this  time  we  can  only  speculate  about  this  and 
outline  its  potential  influence.  Independent  studies  are 
needed  to  develop  a better  understanding  of  this  phenomenon. 

A great  deal  could  be  learned  from  two-dimensional  finite 
difference  calculations  in  which  spallation  is  allowed  to 
occur. 


6.5  Simultaneous  Correction  for  the  Spall  Impulse  and  Double- 
Couple  Contributions  to  the  Source 

! In  Section  V we  discussed  the  double-couple  enhancement 

I 

of  the  source  according  to  the  results  of  Toksoz  and  Kehrer 
(1972).  For  PILEDRIVER  the  double-couple  source  dominates 
' the  Rayleigh  wave,  while  it  plays  a relatively  small  role 

j for  most  of  the  other  events  (Table  7) . In  simultaneously 

j including  all  the  effects  we  have  discussed,  we  must  decide 

what  to  do  with  the  upgoing  waves  generated  by  the  double- 
couple source.  If  this  source  is  spatially  and  temporally 
contiguous  with  the  explosion,  then  both  components  of  the 
source  shoule  be  treated  the  same.  If  not,  we  may  want  to 
l suppress  only  the  upgoing  waves  from  the  explosion  and  leave 

| the  double-couple  unchanged.  Let  us  examine  the  effect  of 

| these  alternatives. 

In  Figure  24  we  show  the  PILEDRIVER  synthetic  seis- 
i mogram  at  TUC  with  a double-couple  of  strength  F = 3.2.  The 

\ same  record  appeared  in  Figure  19.  We  then  suppress  25  percent 

and  50  percent  of  the  upgoing  waves  from  the  total  source.* 

The  effect  is  quite  strong  with  the  amplitude  reduction  being 
proportional  to  the  fraction  of  upgoing  waves  suppressed. 

♦The  formulation  of  the  theory  for  suppressing  the  upgoing 
«eves  is  due  to  D.  G.  Harkrider. 
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The  match  to  the  observed  waveform  at  TUC  is  not 
nearly  as  good  for  PILEDRIVER  as  it  is  for  the  Pahute  Mesa 
and  Yucca  Flat  events,  no  matter  what  we  assume  about  the 
source.  There  seems  to  be  something  about  this  event  that 
is  not  accounted  for  in  our  models. 

If  we  do  not  suppress  any  of  the  upgoing  wave  energy 
from  the  double-couple,  simultaneous  correction  for  double- 
couple and  spallation  effects  is  done  by  a simple  superposition. 
The  double-couple  correction  factors  used  in  Section  V are 
applied  to  the  composite  source  solutions  of  Figures  21  to  23. 
The  results  are  given  in  Table  10. 

Our  final  results  are  summarized  in  Table  10.  We 
first  list  the  mean  source  amplitude  obtained  when  assuming 
the  reduced  displacement  potential  represents  the  entire  source. 
The  mean  double-couple  corrections  deduced  in  Section  V are 
then  listed  as  are  our  preferred  solutions  for  the  spall- 
corrected  source  level.  The  latter  are  from  the  composite 
source  with  25  percent  of  the  upgoing  waves  suppressed  in 
Figures  21  to  23.  The  new  corrected  source  level  is  then 
obtained  by  dividing  the  spall-corrected  source  by  the  double- 
couple correction. 

For  PILEDRIVER  we  give  two  solutions  with  the  first 
obtained  as  outlined  above.  If  we  suppress  the  upgoing  waves 
from  the  double-couple  in  the  same  proportion  as  from  the 
explosion,  the  Rayleigh  waves  are  reduced  as  shown  in  Figure 
24.  The  double-couple  correction  is  reduced  accordingly; 
that  is,  by  two-thirds  for  suppression  of  25  percent  of  the 
upgoing  waves. 

For  the  Yucca  Flat  events,  the  solutions  are  quite 
satisfactory.  For  Pahute  Mesa  events,  the  main  issue  is  the 
extent  of  suppression  of  the  upgoing  waves  by  scattering  and 
spallation.  More  theoretical  work,  particularly  in  the  two- 
dimensional  modeling  of  explosions,  is  required  to  develop 
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TABLE  10 


SIMULTANEOUS  CORRECTION  FOR  THE  EFFECT  OF  SPALLATION 
AND  A SUPERIMPOSED  DOUBLE -COUPLE 


ALQ 

TUC 

f * (TUC) 

(ALQ) 

Yucca  Flat 

Explosion  source  (Table 
Double-couple  correction 

4) 

9.1 

13.2 

1.50 

(Section  V) 

Spall  corrected  source 

1.19 

1.36 

— 

(Figure  21)* 

10.2 

13.9 

1.36 

New  corrected  source 

8.6 

10.2 

1.19 

Pahute  Mesa 

Explosion  source  (Table 
Double-couple  correction 

4) 

6.6 

11.2 

1.57 

(Section  V) 

Spall  corrected  source 

1.62 

1.84 

4> 

(Figure  23)* 

7.8 

14.3 

1.701 

New  corrected  source 

4.8 

7.8 

1.50 

PILEDRIVER 

Explosion  source  (Table 
Double-couple  correction 

4) 

3.9 

6 . 2 

1.59 

(Section  V) 

Spall  corrected  source 

3.55 

3.00 

(Figure  22)  * 

5.0 

7.8 

1.56 

New  corrected  source 

Modified  double-couple  source 

1.4 

2.6 

1.85 

(Figure  23) * 

2.1 

3.8 

1.81 

Based  on  suppression  of  25  percent  of  the  upgoing  waves. 
The  mean  ratio  after  correcting  each  event  individually. 
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a clearer  understanding  of  this  effect.  For  PILEDRIVER,  there 
are  several  aspects  that  remain  puzzling.  The  source  levels 
are  quite  small,  perhaps  too  small  to  be  consistent  with  the 
evidence  from  other  types  of  data  (Table  6) . The  discrepancy 
between  the  ALQ  and  TUC  estimates  is  very  large,  especially 
when  compared  to  that  for  the  other  events.  Finally,  we  have 
the  relatively  poor  match  between  the  synthetic  and  observed 
records  at  TUC.  Considering  all  this,  we  have  much  less 
confidence  in  our  source  estimate  for  PILEDRIVER.  

The  difficulty  with  the  PILEDRIVER  source  estimate 
is  that  it  is  strongly  influenced  by  two  poorly  understood 
effects.  The  explosion  source  level  is  small  because  most 
of  the  Rayleigh  wave  is  from  the  double-couple  if  the  Toksoz 
and  Kehrer  (1972)  solution  is  used.  But  suppression  of  the 
upgoing  waves  has  an  equally  strong  effect,  increasing  the 
inferred  source  level.  At  this  time,  we  are  unable  to  resolve 
these  conflicting  effects,  but  only  point  them  out  and  indicate 
their  importance. 
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VII.  EXAMINATION  OF  THE  DIFFERENCES  BETWEEN 
THE  ALQ  AND  TUC  SOURCE  ESTIMATES 


I 


I 


If  a spherically  symmetric  point  source  is  assumed, 
the  source  level  inferred  from  the  TUC  observations  is  about 
1.5  times  larger  than  that  inferred  from  the  ALQ  observations 
(Table  4) . This  difference  is  quite  consistent;  averaging 
all  22  events  in  the  population  of  Table  4,  the  ratio  is  1.52 
with  a standard  deviation  that  is  25  percent  of  the  mean.  The 
question  is,  why  should  the  source  appear  larger  when  viewed 
from  TUC  than  when  viewed  from  ALQ?  One  possible  explanation 
is  that  there  is  a true  asymmetry  at  the  source  that  causes 
Rayleigh  waves  to  be  more  strongly  excited  toward  TUC.  In 
the  last  two  sections  we  saw  that  correcting  for  the  double- 
couple and  spall  closure  contributions  does  remove  much  of 
the  difference  for  Yucca  Flat.  However,  this  explanation  does 
not  seem  to  work  for  PILEDRIVER  and  Pahute  Mesa. 

A likely  cause  for  differences  in  the  ^ estimates  at 
the  two  stations  is  that  our  plane- layered  models  are  unable 
to  represent  the  laterally  varying  non-planar  earth  with  the 
same  accuracy  for  both  paths.  Lateral  heterogeneities  can 
certainly  cause  amplitude  variations  by  focusing  or  defocus- 
ing  the  surface  waves  at  a particular  site.  Von  Seggern, 
et  al. , (1975)  discuss  focusing  due  to  lateral  variations  in 
Rayleigh  wave  phase  velocity  and  show  that  effects  of  the 
size  seen  here  can  occur.  Focusing  that  would  account  for 
the  asymmetry  in  the  inferred  source  amplitude  must  act  to 
increase  the  observed  amplitudes  at  TUC  and/or  decrease  those 
at  ALQ. 

One  lateral  inhomogeneity  we  know  about  is  the  crustal 
thickening  along  the  NTS-ALQ  path.  It  is  plausible  to  suppose 
that  this  inhomogeneity  causes  the  surface  waves  at  ALQ  to  be 
smaller  than  predicted  by  our  plane-layered  model.  Such  am 
error  would  cause  us  to  underestimate  the  from  the  ALQ  data, 

<30 

leading  to  the  kind  of  asymmetry  we  see. 
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For  our  plane-layered  models  there  are  four  main  factors 
controlling  the  amplitude  of  the  synthetic  seismograms  that 
should  be  examined  more  closely.  There  are  the  dispersion, 
the  choice  of  the  structure  in  the  source  region,  the 
transmission  between  the  source  and  path  structures  and  the 
Q models  for  the  two  paths. 

The  average  dispersion  data  is  matched  about  as  well 
for  the  NTS-TUC  path  as  for  the  NTS-ALQ  path  (Figures  3 and  4)  , 
though  the  dispersion  data  itself  is  better  determined  for 
the  NTS-ALQ  path  (Figure  2) . The  TUC  seismograms  are  dominated 
by  a strong  Airy  phase  at  about  8 seconds,  which  does  make 
this  station  more  difficult  to  analyze.  We  can  see  in  the 
comparison  of  synthetic  and  observed  seismograms  in  Figure  5 
that  subtle  interference  effects  control  the  peak  amplitude 
at  TUC  to  a greater  degree  than  at  ALQ . In  fact,  the  comparison 
suggests  that  our  synthetic  seismograms  probably  underestimate 
the  Airy  phase  amplitude  at  TUC.  However,  while  this  error 
does  contribute  to  the  asymmetry  we  are  discussing,  it  is 
probably  not  much  larger  than  10  percent. 

In  computing  the  theoretical  seismograms,  we  have 
assumed  that  the  amplitude  excitation  is  the  same  for  both 
paths.  Calculations  of  theoretical  seismograms  for  the  two 
paths  differ  only  in  the  dispersion  and  the  transmission 
coefficient,  T ( o>)  . The  source  region  models  were  based  on 
the  structure  for  the  NTS-TUC  path  which  lies  entirely  in  the 
Basin  and  Range  Province,  and  this  is  certainly  more  reasonable 
if  we  want  to  use  the  same  source  region  model  for  both  paths. 
One  could  argue  that  the  NTS-ALQ  synthetics  are  more  accurately 
computed  with  source  region  models  based  on  the  NTS-ALQ  crustal 
structure.  However,  if  we  do  so,  we  increase  the  ALQ  4^ 
estimates  by  only  about  5 percent.  This  is  in  the  right 
direction,  but  is  too  small  to  have  much  significance. 
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The  transmission  coefficient,  T(oi),  is  not  a well 
defined  parameter.  It  is  meant  to  account  for  a gradual 
transition  between  the  source  region  model  and  the  average 
model  for  the  travel  path.  There  are  probably  many  gradual 
transitions  that  are  averaged  in  our  procedure.  As  we  dis- 
cussed in  connection  with  Figure  7,  the  T(u)  are  near  unity 
for  periods  of  primary  interest  and  play  no  significant 
role  in  the  calculation  except  at  short  periods  where  they 
do  help  shape  the  waveform.  This  is,  of  course,  because  our 
source  and  path  structures  differ  greatly  only  near  the  surface. 

Perhaps  the  least  understood  feature  of  the  calculations 
is  the  appropriate  y(m)  or  Q model  to  be  used  for  each  path. 

The  attenuation  models  we  used  (Table  3)  are  related  in  that 
the  Q-8  relationship  is  the  same  along  each  path.  This  leads 
to  y values  that  are  somewhat  larger  for  the  NTS-ALQ  path. 

To  reduce  the  discrepancy  between  the  ALQ  and  TUC  inferred 
'¥oo,  the  difference  between  the  y values  would  have  to  be 
larger  still.  For  example,  if  y(cu)  were  about  10-4  km-1 
larger  at  periods  of  interest  for  the  NTS-ALQ  path  and  the 
same  amount  smaller  for  the  NTS-TUC  path,  the  inferred  average 
TUC  f^/ALQ  ratio  would  drop  from  1.56  to  1.32.  Changes 
in  y of  this  amount  are  plausible,  but  at  this  time  cannot 
be  justified  by  any  independent  information.  The  amplitude 
discrepancy  we  see  may  be  viewed  as  a good  indication  that 
our  y estimate  is  too  low  for  NTS-ALQ  and/or  too  high  for 
NTS-TUC,  but  this  question  cannot  really  be  resolved  without 
more  research. 

In  Sections  V and  VI  we  showed  synthetic  seismograms 
for  a variety  of  source  configurations.  The  match  of  synthetic 
and  observed  seismograms  achieved  for  ALQ  was  excellent.  The 
match  to  the  TUC  observations  was  more  sensitive  to  details 
of  the  source,  but  in  most  cases  was  also  quite  good.  However, 
our  synthetic  records  never  match  the  long  period  (%10  to  12 


seconds)  cycles  just  ahead  of  the  Airy  phase.  Comparison 
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of  theoretical  and  observed  spectra  indicate  that  this  mis- 
match might  be  due  to  there  being  too  little  long  period 
energy  in  the  synthetic  spectra.  To  test  this  hypothesis, 
we  constructed  a y(oj)  which  brings  the  theoretical  and 
observed  amplitude  spectra  into  better  agreement.  The  y(co) 
and  the  resulting  seismogram  for  the  Yucca  Flat  explosion 
source  are  shown  in  Figure  25. 

The  new  seismogram  in  Figure  25  gives  an  excellent 
fit  to  the  observed  record.  If  we  applied  the  same  y(c o)  to  the 
best  composite  source  record  of  Figure  21  (Case  G) , the  fit 
would  be  even  better.  However,  the  peculiar  y ( a>)  , which  is 
peaked  near  9 seconds,  implies  extreme  depth  variations  in  Q 
compared  to  the  models  in  Figure  8 and  may  even  imply  negative 
Q at  some  depths. 

Since  the  fit  to  the  observations  is  enhanced  with 
the  special  y(jj)  in  Figure  25,  we  must  either  suppose  that 
this  y(oj)  is  accounting  for  both  anelastic  and  scattering 
effects,  possibly  with  a frequency-dependent  Q,  or  that  this  y(w) 
is  compensating  for  errors  in  our  specification  of  the  source 
and  path  amplification.  However,  the  kinds  of  models  used 
in  this  study  are  unlikely  to  be  able  to  account  for  spectral 
shaping  like  that  represented  by  the  special  y(co)  . Further, 
the  synthetic  and  observed  spectra  are  in  good  agreement  at 
ALQ,  providing  evidence  that  the  error  is  not  in  the  source 
specification.  We  are,  therefore,  inclined  to  attribute  the 
disagreement  of  the  long  period  portion  of  the  observed  and 
synthetic  records  at  TUC  to  the  inability  of  our  rather  simple 
plane- layered  models  to  perfectly  represent  the  real  earth. 

In  summary,  the  fit  to  the  TUC  seismograms  is  not 
quite  as  good  as  for  the  ALQ  records.  This  is  true  both  for 
the  dispersion  and  for  the  amplitude  spectrum.  We  believe 
the  result  is  to  cause  an  overestimate  of  the  source  amplitude 
required  to  match  the  TUC  data,  though  this  cannot  be  con- 
clusively demonstrated.  The  more  dispersed  seismograms  at 
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ALQ  are  fit  very  well  and  it  is  tempting  to  prefer  the 
inferred  from  these  data.  On  the  other  hand,  the  NTS-TUC 
path  is  in  one  tectonic  province  while  the  known  lateral 
variations  along  the  NTS-ALQ  path  are  much  larger.  These 
lateral  variations  would  probably  cause  us  to  underestimate 
the  source  amplitude  required  to  match  the  ALQ  data.  Finally, 
the  Q models  used  for  the  two  paths  could  easily  include 
errors  that  would  cause  the  amplitude  estimates  to  change  by 
20  to  30  percent.  Changes  much  larger  than  that  seem  unlikely. 

Considering  all  these  factors,  the  difference  in  the 
source  estimates  from  the  two  stations  is  not  surprising.  The 
true  value  is  probably  somewhere  between  the  two,  though  we 
give  somewhat  more  credence  to  the  ALQ  estimates.  As  is  so 
often  the  case,  a better  definition  of  the  answer  can  probably 
only  be  obtained  by  going  through  the  same  process  for  more 
stations  so  random  uncertainties  can  be  reduced  by  averaging. 
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VIII.  CONCLUSIONS 


Our  purpose  in  this  work  was  to  explore  in  as  much 
detail  as  possible  the  usefulness  of  synthetic  seismogram 
techniques  for  studying  underground  nuclear  explosions  and, 
particularly,  for  determining  their  yield.  To  do  so,  we  have 
used  observations  at  two  particular  stations  to  determine 
the  'Pao  for  NTS  explosions.  The  relationship  between  and 
explosion  yield  is  dependent  on  the  source  material  properties 
and  is  a separate  issue  not  discussed  here. 

The  capability  to  determine  the  from  comparison 
of  synthetic  and  observed  seismograms  was  conclusively 
demonstrated.  There  are  uncertainties  and  unresolved  issues, 
but  we  have  attempted  to  carefully  bound  the  range  of  uncer- 
tainty and  to  pose  the  questions  that  remain  unanswered. 

The  key  features  of  this  work  may  be  summarized 
as  follows: 

• Surface  wave  observations  of  explosions  can  be  used 
to  infer  models  for  the  earth  structure  along  the 
travel  path.  These  models  then  predict  synthetic 
seismograms  that  closely  resemble  the  observed 
records . 

• The  amplitude  of  the  source  can  be  inferred  by 
comparing  synthetic  and  observed  records.  For 
Yucca  Flat  and  Pahute  Mesa  events,  the  inferred 
source  amplitude  is  in  general  agreement  with  values 
obtained  by  other  methods . 

• The  effect  of  spallation  was  carefully  explored. 

The  spall  closure  impulse  is  too  small  to  have 
much  influence.  More  likely  to  be  important  is 
the  loss  of  energy  from  the  upgoing  waves  that  is 
associated  with  spallation. 

• The  extent  to  which  the  upgoing  waves  are  attenuated 
by  spallation  or  scattered  is  poorly  understood. 

For  the  wavelengths  of  interest,  it  may  not  even 
occur  to  a signigicant  extent.  More  theoretical 
work  is  required  to  resolve  the  questions. 
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Particularly  useful  would  he  the  study  of  Rayleigh 
waves  from  two-dimensional  finite  difference  cal- 
culations in  which  spall  is  allowed. 


• Tectonic  strain  release  will  affect  the  source 
estimates,  but  it  can  be  corrected  if  sufficient 
azimuthal  coverage  is  available. 

• The  important  PILEDRIVER  event  is,  unfortunately, 
the  most  puzzling  of  those  studied.  The  tectonic 
release  and  spall  related  phenomena  dominate  the 
source  determination.  The  agreement  of  synthetic 
and  observed  waveforms  is  significantly  worse  for 
this  event,  suggesting  that  something  might  be 
unaccounted  for  in  the  vicinity  of  the  source. 

• The  source  level  for  the  Yucca  Flat  events  is  best 
determined.  A substantial  proportion  of  the 
deviation  from  the  mean  for  individual  events 
is  probably  associated  with  real  differences  in 
the  source  coupling.  It  is  important  to  know  how 
large  these  differences  can  be  for  superficially 
identical  events  (individual  ^ values  are  given 
by  Bache,  1978) . 

• The  source  level  for  Pahute  Mesa  events  is  determined 
with  less  confidence  than  for  Yucca  Flat  because 

of  larger  uncertainties  associated  with  tectonic 
release  and  spallation  phenomena.  Also,  relatively 
few  Pahute  Mesa  events  were  examined. 

• Finally,  we  point  out  that  our  data  for  this  study 
was  taken  from  WWSSN  film  clips.  This  is  far  from 
ideal  for  many  reasons,  especially  in  the  narrow 
yield  range  to  which  we  were  restricted.  Digital 

data  from  well-maintained  stations  would  be  much  better. 

Applying  the  techniques  of  this  report  to  foreign 
explosions,  we  would  obtain  a yield  estimate  that  is  independent 
of  the  estimate  from  empirical  Ms-yield  curves.  Of  course, 
the  two  estimates  would  be  the  same  if  the  empirical  data  were 
not  biased  by  some  systematic  factor  that  was  correctly  ac- 
counted for  by  our  more  detailed  modeling  procedure.  The 
uncertainties  in  yield  estimates  from  the  matching  of  synthetic 
and  observed  records  may  be  summarized  as  follows: 
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1.  Different  estimates  would  be  obtained  from  different 
stations.  After  correcting  for  radiation  pattern 
effects,  taking  Love  waves  into  account,  these 
would  be  averaged  and  a statistical  estimate  of 

the  uncertainty  would  result. 

2.  Neither  the  velocity  nor  Q models  will  ever  be 
determined  with  ideal  accuracy.  We  can  estimate 
the  errors  associated  with  uncertainties  in  these 
models  and  the  models  themselves  should  improve 
with  time  if  the  effort  is  made. 

3.  The  most  important  uncertainties  are  associated 
with  the  source  itself.  First,  there  is  the 
relationship  between  the  source  amplitude  and 
the  explosion  yield.  However,  the  most  that  can 
ever  be  determined  with  seismic  methods  is  the 
source  amplitude.  Second,  there  are  the  uncer- 
tainties associated  with  the  selection  of  the 
appropriate  model  for  the  spallation  process  and, 
particularly,  the  loss  of  energy  from  the  upgoing 
waves.  More  research  will  help  resolve  these 
questions.  Also,  we  point  out  that  there  is  no 
reason  to  suppose  foreign  and  U.S.  explosions 
would  be  much  different  as  far  as  this  aspect  of 
the  physics  is  concerned. 

If  we  select  a particular  model  for  this  process 
(for  example,  assuming  that  25  percent  of  the  upgoing  waves 
are  suppressed)  and  interpret  U.S.  and  foreign  events  with 
this  assumption,  our  source  amplitude  estimates  will  be 
consistent. 
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